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Abstract 

The galactic environment of the Sun varies over short timescales as the Sun 
and interstellar clouds travel through space. Small variations in the dynamics, 
ionization, density, and magnetic field strength in the interstellar medium (ISM) 
surrounding the Sun can lead to pronounced changes in the properties of the 
heliosphere. The ISM within ~30 pc consists of a group of cloudlets that flow 
through the local standard of rest with a bulk velocity of 17-19 km s^^, and 
an upwind direction suggesting an origin associated with stellar activity in the 
Scorpius-Centaurus association. The Sun is situated in the leading edge of this 
flow, in a partially ionized warm cloud with a density of ~0.3 cm^"^. Radiative 
transfer models of this tenuous ISM show that the fractional ionization of the 
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ISM, and therefore the boundary conditions of the hehosphere, will change from 
radiative transfer effects alone as the Sun traverses a tenuous interstellar cloud. 
Ionization equilibrium is achieved for a range of ionization levels, depending 
on the ISM parameters. Fractional ionization ranges of x(H)=0. 19-0.35 and 
X(He)=0.32-0.52 are found for tenuous clouds in equilibrium. In addition, both 
temperature and velocity vary between clouds. Cloud densities derived from 
these models permit primitive estimates of the cloud morphology, and the time- 
line for the Sun's passage through interstellar clouds for the past and future ~ 10^ 
years. The most predictable transitions happen when the Sun emerged from the 
near vacuum of the Local Bubble interior and entered the cluster of local inter- 
stellar clouds flowing past the Sun, which occurred sometime within the past 
140,000 years, and again when the Sun entered the local interstellar cloud now 
surrounding and inside of the solar system, which occurred sometime within the 
past 44,000 years, possibly a 1000 years ago. Prior to ~ 140,000 years ago, no 
interstellar neutrals would have entered the solar system, so the pickup ion and 
anomalous cosmic ray populations would have been absent. The tenuous ISM 
within 30 pc is similar to low column density observed globally. In this chap- 
ter, we review the factors important to understanding short-term variations in 
the galactic environment of the Sun. Most ISM within 40 pc is partially ionized 
warm material, but an intriguing possibility is that tiny cold structures may be 
present. 

Keywords: Interstellar Matter, Heliosphere, Equilibrium Models. 

6.1 Overview 

In 1954 Spitzer noted that the "Study of the stars is one of [mankind's] oldest 
intellectual activities. Study of the matter between stars is one of the youngest." 
Comparatively, the study of interstellar matter (ISM) at the heliosphere is an 
infant. Still younger is the study of the interaction between the heliosphere and 
ISM that forms the galactic environment of the Sun. The total pressure of the 
ISM at the solar system is counterbalanced by the solar wind ram pressure, but 
since both the Sun and clouds move through space, this balance is perturbed as 
the Sun passes between clouds with different velocities or physical properties. 
In this chapter, we focus on variations in the galactic environment of the Sun 
over timescales of ±3 Myrs, guided by data and models of local interstellar 
clouds. In essence, we strive to provide a basis for understanding the "galactic 
weather" of the solar system over geologically short timescales, and in the 
process discover recent changes in the Sun's environment that affect particle 
populations inside of the heliosphere, and perhaps the terrestrial climate. The 
heliosphere, the bubble containing the solar wind plasma and magnetic field, 
dances in the wind of interstellar' gas drifting past the Sun. This current of 
tenuous partially ionized low density ISM has a velocity relative to the Sun of 
~26 km s~^. It would have taken less than 50,000 years for this gas to drift 
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from the vicinity of the closest star a Cen, near the upwind direction in the 
local standai^d of rest (LSR, Tables 6.1), and into the solar system. 

Star formation disrupts the ISM. The nuclear ages of massive nearby stars 
in Orion (at a distance of ~ 400 pc) and Scorpius-Centaurus (at ~ 150 pc) 
are 4-15 Myr, so the solar system has been bombarded by high energy photons 
and particles many times over the past 10^ years. Recent nearby supernova 
events include the formation of the Geminga pulsar ~ 250, 000 years ago, 
and possibly the release of ( Oph as a runaway star ~0.5-l Myr ago. The 
dynamically evolving interstellar medium inhibits the precise description of 
the solar Galactic environment on timescales longer than ~ 1 — 3 Myrs. 

Variations in the galactic environment of the Sun for time scales of ~3 Myrs 
are short compared to the vertical oscillations of the Sun through galactic plane, 
and short compared to the disruption of local interstellar clouds by star forma- 
tion. The oscillation of the Sun in the Galactic gravitational potential carries 
the Sun through the galactic plane once every 34 Myrs, and the Sun will 
reach a maximum height above the plane of 78 pc in about 14 Myrs (Bash, 
1986, Vandervoort and Sather, 1993, also see Chapter 5 by Shaviv). The mo- 
tion of the Sun compared to the kinematics of some ensemble of neai^by stai^s 
is known as the solai" apex motion, which defines a hypothetical closed cir- 
cular orbit around the galactic center known as the Local Standard of Rest 
(LSR, §6.1.2). Note that the LSR definition is sensitive to the selection of the 
comparison stars, since the mean motions and dispersion of stellar populations 
depend on the stellar masses. The transformation of the solar trajectory into the 
LSR for comparison with spatially defined objects, such as the Local Bubble, 
introduces uncertainties related to the LSR. 

A striking feature in the solar journey through space is the recent emergence 
of the Sun from a near- vacuum in space, the Local Bubble, with density p < 
10-26 g cm-3 (Frisch, 1981, Frisch and York, 1986, §6.2). We will show that 
the Sun has spent most of the past ^^3 Myrs in this extremely low density region 
of the ISM, known as the Local Bubble. The Local Bubble extends to distances 
of over 200 pc in parts of the third galactic quadrant (/ = 180° — > 270°, 
Frisch and York, 1983), and is part of an extended interarm region between the 
Sagittarius/Carina and the Perseus spiral arms. The proximity of the Sun to 
the Local Bubble has a profound effect on the solar environment, and affects 
the local interstellar radiation field, and dominates the historical solar galactic 
environment of the Sun (§6.2.2). Neutral ISM was absent from the bubble, and 
byproducts of the ISM-solar wind interaction, such as dust, pickup ions and 
anomalous cosmic rays, would have been absent from the heliosphere interior. 

The Local Bubble is defined by its geometry today, so the solar space motion 
is a variable in estimating the departure of the Sun from this bubble. For all 
reasonable solar apex motions, sometime between 1000 and 140,000 years ago 
the Sun encountered material denser than the nearly empty bubble. During the 
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late Quaternary geological period, the Sun found itself in a flow of warm low 
density ISM, n ^ 0.3 cm~^, originating from the direction of the Scorpius- 
Centaurus Association (§§6.3, 6.6.4). The Sun is presently surrounded by this 
warm, T ~ 7, 000 K, tenuous gas, which is similar to the dominant form of 
ISM in the solar neighborhood. The physical characteristics of this cluster of 
local interstellar cloudlets (CLIC) are discussed in §6.3. 

Figure 6.1. SEE figla.pdf and figlb.pdf. The solar position and motion are compared to the 
Local Bubble for two planes aligned with the solar apex motion (Table 6.1). The left figure 
gives the Local Bubble ISM distribution in a plane that is tilted by 27° to the Galactic plane 
with the northern surface normal pointing towards (220°, +63°). Stars with latitudes within 
about ±15° of that plane are plotted. The right figure shows a meridian slice perpendicular 
to the galactic plane, and aligned with the £=40° — 220° axis. Stars with longitudes within 
±25° of the meridian slice are plotted. These ISM distributions were constructed from cleaned 
and averaged photometric and astrometric data for O, B, and A stars in the Hipparcos cata- 
log. For a solar LSR velocity of 13 — 20 pc/Myrs, the Sun has been within the Local Bubble 
for over 3 Myrs. The symbols show E{B — V) values 0.017-0.051 mag (tiny x's), 0.051- 
0.088 mag (boxes), 0.088-0.126 mag (dots), and > 0.126 mag (triangles) so that shading 
levels give A*' (H)~ 10^" '' -> 10^° '', - 10^° '' 10^''■^ and > 10^" ■^ after assuming 
N(H°+2H2 )/-E(-B — V)=21.16 cm~^mag~^. The solar motion is shown for both Standard and 
Hipparcos values. The large intensely shaded triangles show dust clouds from Dutra and Bica 
(2002). The CLIC column densities are below the minimum value of Af(H°)- 10^" "^o cm"^ 
displayed here. 

The solar entry into the CLIC would correspond to the appearance of inter- 
stellar dust, pickup ions and anomalous cosmic rays in the heliosphere. The 
dust filtration by the solar wind and in the heliosheath is sensitive both to dust 
mass and the phase of the solar cycle (see Landgraf chapter), so that these two 
effects are separable. In principle the solar journey through clouds could be 
traced by comparative studies of interstellar dust settling onto the surfaces of 
inner versus outer moons. 

In §6.5.2 we discuss the Sun's entry into, and exit from, the interstellar cloud 
now suiTounding the solai^ system, and show that it can be determined from a 
combination of data and theoretical models. The velocity of this cloud, known 
as the local interstellar cloud (LIC), has been determined by observations of 
interstellar He° inside of the solar system (Witte, 2004). We test several possi- 
bilities for the three-dimensional (3D) shape of the LIC, and conclude that the 
Sun entered the surrounding cloud sometime within the past 40,000 years, and 
possibly very recently, in the past '^lOOO years. 

The velocity difference between interstellar He° inside of the solar system, 
and interstellar D°, Fe"*" and other ions towards stars in the upwind direction, 
including the nearest star a Cen, together indicate the Sun will emerge from the 
cloud now surrounding the Sun within ~4,000 years. This raises the question 
of "What is an interstellar cloud?", which is difficult to answer for low column 
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density ISM such as the CLIC, where cloudlet velocities may blend at the 
nominal UV resolution of 3 km s^^. 

Determining the times of the solar encounter with the LIC, and other clouds, 
requires a value for n(H°) in the cloud because distance scales are determined 
by A^(H°)/n(H°). Theoretical radiative transfer models of cloud ionization pro- 
vide n(H°) and n(H"'") for tenuous ISM such as the LIC, and supply important 
boundary constraints for heliosphere models. CLIC column densities are low, 
N ^ 10^^, and radiative transfer effects are significant (§6.4). Models pre- 
dict n(H°)~0.2 cm"'^ and n(H~'^)~0. 1 cm"'^ at the Sun, with ionization levels 
rising slowly to the cloud surface so that (n(H°))Lic~0-17 cm^'^. 

In §6.6.1 we discuss the relative significance of diffuse warm partially ion- 
ized gas (WPIM) and traditional Stromgren spheres for the heliosphere. Dis- 
tinct regions of dense fully ionized ISM are infrequent near the Sun, although 
regions of ionized gas surround distant massive stars bordering the Local Bub- 
ble, and surround (3 Cen and A Sco, which helps energize the Loop I superbub- 
ble interior. The Sun's path is unlikely to traverse, or have traversed, an H II 
region around an 0-B 1 star for ±4 Myrs because no such stars or H II regions 
are within 80 pc of the Sun. Hot white dwarf stars are more frequent nearby, 
and may in some cases be suiTounded by small Stromgren spheres. In con- 
trast, diffuse ionized gas is widespread and has been observed extensively us- 
ing weak optical recombination lines from H"*" and other species, synchrotron 
emission, and by the interaction of pulsai^ wave packets with ionized interstel- 
lai- gas. The CLIC and diffuse ionized gas are part of a continuum of ionization 
levels predicted and observed for low density interstellar gas. 

The heliosphere should vary over geologically short timescales, for reasons 
unrelated to the Sun or solar activity, when encountering the small scale struc- 
ture of the nearby ISM. These variations in the cosmic environment of the 
Sun may be traced by terrestrial radioisotope records, or by variable fluxes 
of interstellar dust grains deposited on planetary moons (see other chapters in 
this volume). Over 3 Myrs timescales the Sun travels 40-60 pc through the 
LSR, and the interval between cloud traversals should be less than < 80, 000 
years. Differing physical properties of cloudlets in the flow will affect the 
heliosphere and interplanetary medium. If cloud densities are relatively uni- 
form at ri(H°)~0.2 cm~'^, then data on ISM towards nearby stars indicate these 
clouds fill ~33% of space within 10 pc, with mean cloud lengths of 1.0 pc, 
and with a downwind hemisphere that is emptier than the upwind hemisphere 
(§6.3.3). 

The ISM within ~2 pc is inhomogeneous. The ratio Fe+/D° increases as 
the viewpoint sweeps from the downwind to upwind direction, signaling either 
grain destruction in an upwind interstellar- shock, or increased cloud ionization. 
The Fe"'"/D° gradient is consistent with the destruction of CLIC dust grains 
in a fragment of the expanding superbubble around the Scorpius-Centaurus 
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Association (§6.6.4). There are several possibilities for the next cloud to be 
encountered, and the most likely candidate is the cloud towards the nearest 
star Q Cen (§6.5.3), where only modest consequences for the heliosphere are 
expected, since the temperature and velocity of the cloud differ only by ~800 
K and ~3 km s~^, respectively. 

The properties of the CLIC are similar to low column density ISM observed 
in the solar vicinity, and observed in the Arecibo Millennium survey of very 
low A^(H°) clouds (Heiles and Troland, 2004). The warm neutral medium 
(WNM) in the Arecibo survey appears to be similar to the CLIC, but so far, 
counterparts of the exotic cold, low column density, neutral clouds, T ~ 50 
K and A^(H°)~ 10^^ cm^^, are not found locally (§6.6). If they are hidden in 
the upwind ISM, they would be far from ionization equilibrium (having much 
lower ionization than the surrounding warm, ionized gas), in contrast to the 
Lie gas (§6.4). Cold ISM towards 23 Ori indicate densities of >10 cm~^ for 
the cold neutral medium gas (CNM). 

We do not discuss high-velocity, v ~ 100 km s~^, low column density ra- 
diative shocks. Such disturbances would cross the heliosphere on timescales 
of 5 — 10 years, seriously perturbing it (Sonett et al., 1987, Frisch, 1999, 
Mueller et al., 2005), and would be evolving as they cool. Low column density 
high velocity shocked gas, A^(H)~ 10^^ — 10^^ cm~^, no~0.5-5 cm~^, could 
be common in ISM voids. However, with high velocities of ~100 km s~^, 
such gas will traverse voids such as the Local Bubble in less than 1 Myr. Sim- 
ilar high-velocity gas is associated with the Orion's Cloak superbubble shell 
around Orion and Eridanus, with v ~ 100 km s"^ and 7V(H)~ 10^^ - 10^^ 
cm"^ (Welty et al., 1999). Rapidly moving low column density shocks similar 
to these would pass over the heliosphere quickly, in 50 - 1,500 years, with 
strong impUcations for heliosphere structure (Zank et al., this volume). The 
fairly coherent velocity and density structure and relatively weak turbulence in 
the CLIC argues against any such disturbance having occurred recently (see 
§6.3.4). We also do not discuss dense molecular clouds, which are discussed 
in the Shaviv chapter. 

As a way of summarizing short term solar environment, we show the solar 
position inside of the Local Bubble in Figure 6.1 (see 6.2.1). We will now 
briefly review relevant fundamental material. 

6.1.1 Fundamental Concepts 

Basic concepts needed for understanding the following discussions are sum- 
marized here. We use the local standard of rest, LSR, based on the Standard so- 
lar motion, unless otherwise indicated. Velocities quoted as heliocentric (HC) 
are with respect to the Sun, while LSR velocities are obtained by correcting 
for the motion of the Sun with respect to the nearest stars, or the solar apex 
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motion. The basic principles underiying obtaining the properties of interstellar 
clouds are summarized in §6.1.3. 

The ultraviolet and X-ray parts of the spectrum play an essential role in 
the equilibrium of interstellar clouds. The spectral regions relevant to the fol- 
lowing discussions are the ultraviolet (UV), with wavelengths 1200 A < A < 
3000 A; the far ultraviolet (FUV), 912 A < A < 1200 A; the extreme ultravio- 
let (EUV), 100 A < A < 912 A; and the soft X-ray, A < 100 A. The spectrum 
of the soft X-ray background (SXRB), photon energies < 0.5 keV, has been 
an important diagnostic tools for the distribution of ISM with column density 
A'^(H°)> 10^^'^ cm~^. The first ionization potential (FIP) of an atom is the 
energy required to ionize the neutral atom. The FIP of H is 13.6 eV. Elements 
with FIP < 13.6 eV are predominantly ionized in all clouds, since the ISM 
opacity to photons with energies E < 13.6 eV is very low. The ionization 
potential, IP, of several ions are also important. For instance IP(Ca+)=l 1.9 eV, 
and iV(Ca++)/Af(H°) may be large in warm clouds. For T > 4, 000 K and 
n(e)<0.13 cm~'^, Ca+^/Ca"*" > 1 . Variations in the spectrum of the ionizing 
radiation field, as it propagates through a cloud caused by the wavelength de- 
pendence of the opacity, r(A), are referred to as radiative transfer (RT) effects. 

Color excess is given by E{B — V), and represents the observed star color, 
in magnitudes, measured in the UBV system, compared to intrinsic stellai^ col- 
ors (e.g. Cox, 2000). Color excess represents the intrinsic star color combined 
with reddening by foreground dust. The column density, A^(X), is the number 
of atoms X contained in a column of material with a cross section of 1 cm~^ 
and bounded at the outer limit by an arbitrary point such as a star. The volume 
density, n{X), is the number of atoms cm""^. The notation nHi,o.2 indicates the 
H° space density in units of 0.2 cm~^. The fractional ionization of an element 
X is N{X+)/iN{X°)+N{X+)). The filling factor, /, gives the percentage of 
the space that is filled by ISM denser than the hot gas in the Local Bubble (i.e. 
warm or cold gas). 

Among the various acronyms we use are the following: The interstellar 
cloud feeding neutral ISM and interstellar dust grains (ISDG) into the solar 
system is the local interstellar cloud (LIC). A cluster of local interstellar clouds 
(CLIC) streams past the Sun, forming absorption lines in the spectra of neai^by 
stai^s. Alternative names for the CLIC ai^e the "Local Fluff" and the very local 
ISM, and the LIC is a member of the CLIC. The upwind portion of the CLIC 
has been called the "squall line" (Frisch, 1995). In the global ISM we find 
cold neutral medium (CNM), warm neutral medium (WNM), and warm ion- 
ized medium (WIM). Both observations and our radiative transfer models of 
tenuous ISM, A^(H) ^ 10^^ cm~^, indicate that warm partially ionized medium 
(WPIM) is prevalent nearby. 

The Scorpius-Centaurus Association (SCA) is a region of star formation 
that dominates the distribution of the ISM in the upwind direction of the CLIC 
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flow, including the Loop I supernova remnant. The Local Bubble (LB) is a 
cavity in the interstellar dust and gas around the Sun that may merge in some 
regions with the bubble formed by stellar evolution in the SCA, Loop L near 
the solar location (Frisch, 1995). The Local Bubble is used here to apply to 
those portions of nearby space, within ~70 pc, with densities p < 10~^^ g 
cm~^. 

The densities of neutral interstellar gas at the solar location are provided by 
data on pickup ions (PUI) and anomalous cosmic rays (ACR). PUIs are formed 
from interstellar neutrals inside the heliosphere, which have become ionized 
and captured by the solar wind (Moebius et al., this volume). ACRs are PUIs 
that are accelerated in the solar wind and termination shock to energies of < 1 
GeV/nucleon (Fahr et al., Florinski et al., this volume). 

Element abundances in the CLIC gas are similar to warm gas in the Galac- 
tic disk, and the missing atoms are believed to be incorporated in dust grains 
(Frisch et al., 1999, Welty et al., 1999). We will assume that the solar abun- 
dance pattern provides the benchmark abundance standard for the local ISM, 
although this issue has been debated (Snow, 2000), and although the solar 
abundance of O and other elements are controversial (e.g. Lodders, 2003). 
Using solar abundances as a baseline, the depletion 5 of an element repre- 
sents the amount of the element that is missing from the gas. Depletion is 
defined as 5x = log [xig^f" element X with abundance [X] (Savage 

and Sembach, 1996). The depletion and condensation temperature of an el- 
ement correlate fairly well. This is shown by the three observed interstellar 

depletion groups 5Ti,Ca ' 3.3, 5Fe,Cr,Co,Niy 2.3, and 6Mg,Si ~ -1-3, 

which condense at temperatures of ~ 1, 500 K, 1, 330 K, and ^ 1, 300 K 
respectively during mineral formation in an atmosphere of solar composition 
and pressure (Ebel, 2000). The interstellar gas-to-dust mass ratio and grain 
composition can be reconstructed from the atoms missing from the gas phase 
(Frisch et al., 1999). Depletions vary between warm tenuous and cold dense 
ISM, an effect originally parameterized by velocity and known as the Routly- 
Spitzer effect (McRae Routly and Spitzer, 1952), and now attributed to par- 
tially ionized warm ISM that has been shocked (§6.4.3). The largest variations 
are found for refractory elements, and result from ISDG destruction by inter- 
stellar shocks (Jones et al., 1994, Slavin et al., 2004). Depletion estimates are 
sensitive to uncertainties in A'^(H), which in cold clouds must include A^(H2) , 
and in warm tenuous clouds must include H"^, both of which can be hard to 
directly measure. 

Satellites that have revolutionized our understanding of the galactic envi- 
ronment of the Sun include Copernicus, the International Ultraviolet Explorer, 
lUE, the Extreme Ultraviolet Explorer, EUVE, the shuttle launched Interstel- 
lar Medium Absorption Profile Spectrograph, IMAPS, the Far Ultraviolet Ex- 
plorer, FUSE, and the Hubble Space Telescope, HST (with the Goddard High 
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Resolution Spectrometer, GHRS, and the Space Telescope Imaging Spectro- 
graph, STIS). 

Several papers in the literature have provided important insights to the phys- 
ical properties of the local ISM, and are referred to throughout. Among them 
are: Frisch (1995, hereafter F95), which reviews the physical characteristics 
of the Local Bubble and the CLIC; Gry and Jenkins (2001, hereafter GJ) and 
Hebrard et al. (1999), which present the physical characteristics of the two 
nearby interstellar clouds observed in the downwind direction, the LIC and 
blue-shifted cloud; Slavin and Frisch (2002, hereafter SF02) and Frisch and 
Slavin (2003, 2004, 2005, hereafter FS), which present detailed radiative trans- 
fer model calculations of the surrounding ISM; Frisch et al. (2002, hereafter 
FGW), which probes the kinematics of nearby ISM; and a series of papers by 
Redfield and Linsky (2000, 2002, 2003, 2004, hereafter RL,), and Wood et al. 
(2005, hereafter W05), which assemble a wide range of UV data on the CLIC. 
Recent papers presenting results from FUSE also contribute significantly to 
our understanding of the tenuous ISM close to the Sun. 

6.1.2 The Solar Apex Motion and Local Standard of Rest 

Comparisons between the solar position and spatially defined objects, such 
as the LB, require adoption of a velocity frame of reference. The "Local Stan- 
dard of Rest" (LSR) represents an instantaneous inertial frame for a corotating 
group of nearby stars (<500 pc) on a closed circular orbit around the galactic 
center, and is commonly used for this purpose. The solar motion with respect 
to the LSR, known as the solar "apex motion", is dynamically defined with 
respect to a sample of nearby stars. All reasonable values for the LSR indicate 
that the Sun is traveling away from the ISM void in the third quadrant of the 
Galaxy, £=180°^ 270°, at a velocity of ~ 13-20 pc/Myrs. 

The Sun is ~8 kpc from the galactic center. The mean orbital motion around 
the galactic center of nearby stars is ~220 km s~^, and the solar velocity is 
225ib20 km s^^ (e.g., see proper motion studies of extragalactic radio sources 
beyond the galactic center, Reid et al., 1999). Accurate astrometric data for 
stars is provided by the Hipparcos catalog of proper motions and positions of 
~ 120,000 nearby stars (Ferryman, 1997). For the coordinate system x, y, z, 
which represent unit vectors towards the galactic center, direction of galactic 
rotation, and north galactic pole respectively, the speeds V^, Vy, and Vz then 
represent the mean solar velocity in the x, y, and z, directions compared to 
some arbitrarily selected set of nearby stars. A kinematically unbiased set 
of stars from the Hipparcos Catalog shows that the mean vertical and radial 
components of nearby star velocities have no systematic dependence on stellar 
mass (Dehnen and Binney, 1998). However Vy determined for stars hotter than 
AO {B — V <0.0 mag) exceeds the mean LSR by ~6 km s~^, and the velocity 
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dispersion for these relatively young stars is smaller than for evolved low mass 
stai^s. Since molecular clouds may share the motions of massive stars, the 
result is a lack of clarity concerning the best LSR transformation for studying 
the solar motion with respect to spatially defined objects. 

The general practice is to transform radio data into the LSR velocity frame 
using the "Standard" solar motion, which is based on the weighted mean veloc- 
ity for different populations of bright nearby stars irrespective of spectral class. 
The velocity of the Sun with respect to the Standard LSR (LSR^f^) is {V^, Vy, 
V;)=(10.4, 14.8, 7.3) km s^^ (Mihalas and Binney, 1981), giving a solar apex 
motion, of 19.5 km s~^ towards /=56°, b=+2T (or LSR^td, Table 6.1). The 
solar motion with respect to a kinematically unbiased subsample of the Hip- 
parcos catalog (omitting stars with B — V < 0.0 mag) is (T4,yy,yz)=( 10.0,5.3, 
7.2) km s~^, corresponding to a solar speed of V=13A km s~^ towards the 
apex direction (.=-21. 1° and 6=32.4° km s~^ (or LSR^ip)- 



Table 6.1. The LSR Velocities of the Sun and Nearby Interstellar Clouds 





HC 


LSR^td 


LSRj/jp 




V,ib 


V,£b 


V,£b 




(km s-i, ° °) 


(kms-\ ° °) 


(kms-\ ° °) 


Sun 




19.5, 56° 23° 


13.4, 27.7° 32.4° 


CLIC 


-28.1±4.6, 12.4° 11.6° 


-19.4, 331.0° -5.1° 


-17.0, 357.8° -5.1° 


Lie 


-26.3±0.4, 3.3° 15.9° 


-20.7, 317.8° -0.5° 


-15.7, 346.0° 0.1° 


GC 


-29.1,5.3° 19.6° 


-21.7, 323.6° 6.3° 


-17.7, 351.2° 8.5° 


Apex 


-35.1, 12.7° 14.6° 


-24.5, 341.3° 3.4° 


-23.3,5.5° 4.1° 



Notes: The first row lists the Standard and Hipparcos values for the solar 
apex motion (§6.1.2). The LIC, CLIC, GC (G-cloud), and Apex velocity 
vectors refer to the upwind direction. Columns 2, 3, and 4 give the heliocentric 
velocity vector, and the LSR velocities using coiTcctions based on the Standard 
and Hipparcos-derived solar apex motions, respectively. 



6.1.3 Finding Cloud Physics from Interstellar Absorption 
Lines 

Sharp optical and UV absorption lines, formed in interstellar clouds between 
the Sun and nearby stars, are the primary means for determining the physical 
properties of the CLIC. UV data provide the best look at cloud physics because 
resonant absorption transitions for the primary ionization state of many abun- 
dant elements fall in the UV portion of the spectrum. However, results based on 
UV data are limited by a spectral resolution, which typically is <L 3.0 km s~^, 
and are unable to resolve the detailed cloud velocity structure seen in high res- 
olution (typically ~0.3-0.5 km s~^) optical data. Ultimately high resolution 
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UV data, R ~300,000, are required to probe the CLIC velocity structure and to 
separate the thermal and turbulent contributions to absorption line broadening. 
The details of interpreting absorption line data can be found in Spitzer (1978); 
see alternatively the classic application of the "curve of growth" technique to 
the ISM towards C Oph (Morton, 1975). 

The classic target objects for ISM studies are hot rapidly rotating O, B and 
A stars, where sharp interstellar absorption features stand out against broad 
stellar lines. However O and B stars are relatively infrequent nearby, and in- 
terstellar line data can be contaminated by sharp lines formed in circumstellar 
shells and disks around some A and B stars. Hence care is required to separate 
ISM and circumstellar features (e.g. Ferlet et al., 1993). White dwarf stars of 
DA and DO spectral types provide a hot far UV continuum that can, in many 
cases, be observed in the interval 912-1200 A (e.g. Lehner et al., 2003). Cool 
stars are very frequent, and there are ~ 30 times more G and K stars than A 
stars. The cool stars have a weak UV flux and active chromospheres. Strong 
chromospheric emission lines can be used as continuum sources for interstel- 
lai- absorption lines, although great care is required in the analysis (e.g. Wood 
et al., 2005). The disadvantage of cool stars is that uncertainties arise from 
blending of the interstellar absorption and stellar chr^omosphere emission fea- 
tures. For example, only the Sun has an unattenuated H° Lyman-a emission 
feature, and its strength and shape vary with the phase of solar magnetic activ- 
ity cycle. Heeding these limitations, cool stars provide an important source of 
data on local ISM. 

Absorption Line Data. Observations of interstellar absorption lines in 
the interval 912-3000 A are required to diagnosis the physical properties of 
the ISM, such as composition, ionization, temperature, density, and depletions 
(e.g., Jenkins, 1987, York, 1976, Snow, 2000, Savage, 1995). Most elements 
in diffuse interstellar- clouds ai^e in the lowest electronic energy states, and 
have resonant absorption lines in the UV and FUV The type of ISM traced 
by H°, D°, H2, C+, C+*, N°, N+, 0°, 0+, Mg°, Mg+, Si+, Si+^, Ar°, 
S~^, Fe+ depends on the element FIP and interaction cross-sections. Neutral 
gas is ti-aced by D°, H°, N°, 0°, and Mg°, although Mg° is also formed by 
dielectronic recombination in WPIM (§6.3.1). Charge exchange couples the 
ionization levels of N, O, and H. The CLOUDY code documentation gives 
sources for chai^ge exchange and other reaction rate constants (§6.4). Ions with 
FIP<13.6 eV, such as C+, Mg+, Si+, S+, and Fe+, are formed in both neu- 
tral and ionized gas. Highly ionized gas is traced by Si"^^ and C^^. The ra- 
tios iV(Mg+)/iV(Mg°) and iV(C+)/iV(C+*) are valuable ionization diagnostics 
(§6.3.1). 

Local ISM has been studied at resolutions of 2.5 — 20 km s~^ for ~ 30 
years by Copernicus, lUE, HST, HUT, EUVE, IMAPS, FUSE, and various 
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rocket and balloon experiments. The best source of UV data in the 1190-3000 
A region are data from the GHRS and STIS on the Hubble Space Telescope, 
while FUV data are available from FUSE, IMAPS, and Copernicus. 

Ground-based optical observations have an advantage over space data in two 
respects — instruments provide very high resolution, ~0.3-1.0 km s~^, and 
there are fewer constraints imposed on telescope time. Lines from Ca"*", Na° 
and Ti+ are useful diagnostics of the nearby ISM. The Ca 11 K line (A3933 A) 
is usually the strongest interstellar feature for nearby stars because Ca deple- 
tion is lower (abundances are higher) in warm gas (e.g. Savage and Sembach, 
1996), and the lines are strong. Calcium and titanium are refractory elements 
with highly variable abundances. Both Ca^" and Na° are trace ionization states 
in the local ISM, while Ti+ (F1P=13.6 eV) and H° have similar ionizations. 
High resolution optical data provide the best diagnostic of the cloud kinemat- 
ics of ISM near the Sun and in the upwind direction. Optical lines, including 
Call, are generally too weak for observations in the downwind direction. 

Interpreting Absorption Data. Interstellar absorption features in a stellar 
spectrum provide a fundamental diagnostic of ISM physics. Cloud tempera- 
tures ai^e determined from these data, but the ~3 km s~^ resolution of the best 
UV data is inadequate to separate out closely spaced velocity components, or 
fully distinguish turbulent and thermal line broadening. 

The wavelength-dependent opacity of an absorption line traces the atomic 
velocity distribution and the underlying strength of the transition between the 
energy levels, the Einstein transition probability. The atomic velocity distribu- 
tion has contributions from the bulk ISM motion, turbulence, and the thermal 
dispersion of atomic velocity. In principle, thermal and turbulent contributions 
to the line broadening can be separated, given data for atoms with a range of 
atomic mass. In practice, data on either D° or H°, as well as heavier elements, 
ai^e needed. The temperature of the cloud is found from line widths, so under- 
standing the limitations of this method ai^e important. 

Absorption line data ai^e analyzed with the assumption that particles, of ve- 
locity V and mass ttia, obey a Maxwellian distribution at a kinetic temperature 

f{v) ~ Tj^'^^^ exp(— mAf^/2A;TK). The kinetic temperature also param- 
eterizes the energy of ion-neutral and neutral-neutral elastic scattering. The 
velocity dispersion is characterized by the Doppler line broadening parame- 
ter, which is 6d= (^~")^^^ for ^ purely thermal distribution of velocities. 
The line full-width-at-half-maximum is FWHM~1.7 6d- Few interstellar- ab- 
sorption lines ai^e Maxwellian shaped, however, because the line opacity is a 
non-linear function of column density, and because multiple unresolved clouds 
may contribute to the absorption. In this case, lines are interpreted using a best- 
fitting set of several velocity components, each representing a sepai^ate cloud 
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or "velocity component", at velocity v, that are found to best reproduce the line 
shape given the instrumental resolution. 

Each velocity component, however, is broadened by both thermal motions 
and turbulence. The non-thermal line broadening is incorporated into a general 
term "turbulence", ^, so that = (2A;TK/ma) + The turbulence term 
^ is supposed to trace the mass-independent component of line broadening. 
The main technique for separating turbulence from thermal broadening uses 
6d as a function of atomic mass, although this technique does not distinguish 
unresolved velocity structure. 

For low opacity lines, column density is directly proportional to the equiva- 
lent width W , and N ~ W/\^. For high line center opacities, Tq oc NX/bi), 
the line strength increases in a non-linear way, W ~ A6D(lnro)^/^. Column 
densities derived from observations of partially saturated lines are highly in- 
sensitive to line strength. These important limitations in determining the com- 
ponent properties from absorption line data are well known (Spitzer, 1978). 

Resolution Limitations: Unresolved Clouds. The resolutions of UV 
instruments, R < 10^, are well below those of the best optical spectrometers, 
i? > 3 X 10^. The result is the loss of information about the velocity structure of 
the clouds forming UV absorption lines. High resolution, 0.3-0.5 km s~^, op- 
tical data on Ca+, Na°, and K° show that the number of adjacent components 
separated by velocity 6v increases exponentially as 5v —>■ km s~^ (Welty 
et al., 1994, Welty et al., 1996, Welty and Hobbs, 2001). The result is that 
~60% of cold clouds may be missed at the ~3 km s~^ resolution of STIS and 
GHRS. Cold clouds have a median 6D(Na°)'^0.73 km s~^, and typical temper- 
ature ~80 K. Fig. 6.2 shows the distribution of component separations, Sv, for 
Na°, K°, and Ca"*". The reduction of the numbers of components with small 
separations, 6v ^ 1.5 km s~^, compared to the best fit exponential distribution, 
indicate unresolved velocity structure, while the increase of component sepa- 
rations for 6v > 6 — 7 km s~^ indicates that the distribution of cloud velocities 
is not purely Gaussian (Welty and Hobbs, 2001) 

6.2 Solar Journey through Space: The Past 10^ to 10« 
Years 

The Sun has spent most of Quaternary era, the past 2-3 Myrs, in the Local 
Bubble, a region of space with very low ISM density, which extends to dis- 
tances of over 200 pc in parts of the third galactic quadrant (£=180°— >270°, 
Frisch and York, 1986). Sometime within the past ~ 140,000 years the Sun 
encountered a substantially denser region, and the galactic environment of the 
Sun changed dramatically. The Local Bubble appears to be part of an ex- 
tended interarm region between the Sagittarius/Carina and the Perseus spiral 
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Figure 6.2. The distribution of velocity separations (5v) between adjacent absorption com- 
ponents for Na°, K°, and Ca"*". The dotted line shows the best fit over the range 2 kms"^< 
5v < 6.5 km s~^, with slopes -0.9, -0.6, and -0.7, for K°, Na°, and Ca^ respectively. The 
turnover in the numbers of components for small separations is the result of unresolved velocity 
structure, while the increase of component separations for 5v > 6 — 7 km s~^ shows that the 
distribution of cloud velocities is not purely Gaussian (Welty and Hobbs, 2001). Figure courtesy 
of Dan Welty. 



arms (Beck, 2001). For all viable solar apex motions, the Sun is leaving the 
emptiest part of this bubble and now entering low density ISM approaching us 
from the direction of the Scorpius-Centaurus Association (SCA, Frisch, 1981, 
1995). The solar motion with respect to the LB and to the CLIC are shown in 
Figs. 6.1 and 6.5. The LB dominates the historical solar galactic environment 
of the Sun and affects the local interstellar radiation field. The absence of in- 
terstellar neutrals in the near void of the LB interior would lead to an absence 
of pickup ions and anomalous cosmic rays in the heliosphere for periods when 
the Sun was in the bubble. Reconstructing the solar galactic environment for 
the past several million years requires knowledge of both solar motion and the 
LB properties. 



6.2.1 



Inside the Local Bubble 



The LB surrounding the Sun was discovered through the absence of in- 
terstellar dust grains, which redden starlight, in the nearest ~70 pc (Eggen, 
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1963, Fitzgerald, 1968, Lucke, 1978, Vergely et al, 1997). Reddening data is 
typically sensitive to color excesses E{B — y)>0.01 mag, corresponding to 
A^(H)= 10^^ ''^ cm~^, depending on photometric quality. The LB overlaps the 
interior of Gould's Belt (Frogel and Stothers, 1977, Grenier, 2004). About 1% 
of the ISM mass is in dust grains, and interstellar gas and dust are well mixed 
in space. 

Data showing starlight reddening by interstellar dust are used here to deter- 
mine the Local Bubble configuration (LB, Figs. 6. 1), Fig. 6. 1 displays this con- 
figuration for two planes aligned with the solar apex motion (Table 6.1). The 
right figure shows a slice perpendicular to the galactic plane, where only stars 
with longitudes within 25° of the plane are plotted. The left figure shows the 
distribution of observed warm gas in a plane that is tilted by 27° to the Galactic 
plane with the northern surface normal pointing towards {l,b) = (220°, +63°) 
(thus intersecting the galactic plane along ^=130° and ^=310°). Only stars 
with latitudes within about ibl5°of that plane are plotted. The E{B — V) val- 
ues are smoothed for stai^s with overlapping distances and within 13° of each 
other. Fig. 6.1 shows that in no direction is a column density of A^(H)> 10^^'^ 
cm~^ identified close to the Sun in the anti-apex direction (for the gas-to-dust 
ratio N{Yi)/E{B -V) = 5.8 x lO^^ atoms cm-^ mag-\ Bohlin et al., 1978). 
The reddening data are from color excess values, E{B — V), calculated for O, 
B, and A0-A3 stars in the Hipparcos catalog, after cleaning the data to omit 
variable or otherwise peculiar stars (Hipparcos flag H6=0, Penyman, 1997). 
Both projections show that an ISM column density of A^(H)> 10^^'^ cm~^ is 
not reached in the nearest ^■60-70 pc of the Sun along the solar trajectory, and 
this conclusion is unaltered by the selection of the solar apex motion. 

The Local Bubble cavity is found in optical, UV, and EUV data, with differ- 
ent sensitivities to the ISM. The LB is seen in UV data sensitive to A^(H)> 10^ 
cm~^ (Frisch and York, 1983), in the distribution of white dwarf and cool stars 
observed in the EUV and sensitive to A^(H)> 10^^ cm^^ (Warwick et al., 
1993), in Na+ absorption line data sensitive to A^(H) ^ 10^^ cm^^ (Sfeir 
et al., 1999, Vergely et al., 2001), and in dust polarization data sensitive to 
A^(H) ^ 10^^ cm~^ (Leroy, 1999). The sensitivities and sample densities dif- 
fer between surveys, but the results are consistent in showing that there is no 
slowly moving dense ISM in the path of the Sun for the previous 3 Myrs, or 
~60 pc. UV and EUV data are more sensitive to low column densities, and 
detect upwind CLIC gas. Figure 1 in Chapter 2 displays the Local Bubble 
compared to nearby stellar associations, the Gum Nebula, which is a giant H II 
region, and the distribution of molecular clouds in the solai^ neighborhood. 

Are there nearby interstellar clouds that are inside of the Local Bubble but 
not shown in Fig. 6.1? CLIC column densities are too low for detection in 
E{B — V) data, but the CLIC is seen in EUV, UV, and optical data as local 
ISM with A^(H) ^ 10^^ cm~^ (§6.3). Looking towards the third galactic quad- 
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rant, the only ISM close to the anti-apex direction, and in the distance interval 
5< d <80 pc, are Ca+ components in 7 Ori (^=197°, 6=-16°, d=75 pc) show- 
ing iV(H°)~ lo^^-so cm-2 if Ca+/H°~10~^. The LSR motions of these Ca+ 
components correspond to 0— >13 km s~^(Standard LSR), so that the Sun has 
moved 20-33 pc with respect to this ISM over ^^3 Myrs. Only if one of the 
components towards 7 Ori is within 35 pc of the Sun will they have surrounded 
the Sun recently. This argument is consistent with Ca"*" limits towards the star 
j3 Eri (27 pc), near the anti-apex direction, indicating A^(H°)< lO^^-^*^ cm~^ 
(Frisch et al., 1990). There is not enough data to rule out encounters with very 
low column density clouds near the anti-apex direction, iV(H°) ^ 10^^ cm~^, 
particularly since the Arecibo survey shows that such ISM may have high LSR 
velocities (§6.6). 

Data on CO molecular clouds indicate the Sun is unlikely to have crossed 
paths with a large dusty cloud over the past 3 Myrs. Nearby CO clouds, d < 130 
pc, typically have low LSR velocities. Several dust clouds ai^e found at ~120 
pc (lai^ge triangles. Fig. 6.1) towards the anti-apex direction, £~200°, fe~-36°. 
They ai'e 3C 105.0 (y=8.0 km s"!), MBM20 (^=0.3 km s'^), and LDN1642 
(y=1.3 km s~^. Standard LSR motions). The clouds would have traveled less 
than 24 pc over 3 Myrs, compared to the distance traveled by the Sun of 40 — 60 
pc, so Sun-cloud separations remained over 40 pc. 

6.2.2 Radiation and Plasma of the Local Bubble 

The location of the Sun inside of a cavity in the ISM has a profound affect 
on the immediate solar environment, including the ambient radiation field and 
the properties of the interstellar plasma at the solar location. Prior to the entry 
of the Sun into the CLIC, the very low density material of the Local Bubble 
surrounded the Sun. We now go into some detail about the physical charac- 
teristics of that material, and the resultant radiation field, because the Local 
Bubble radiation field is the dominant factor in the ionization state of the ISM 
close to the Sun, and because the physical properties of the heliosphere are 
extremely sensitive to the LIC ionization. The following discussion is based 
on equilibrium assumptions for the Local Bubble plasma. Alternate models, 
including models of non-equilibrium cooling, are found in the literature but 
not discussed here. 

While the nature and origin of the Local Bubble is the subject of ongoing de- 
bate, at least a few facts are agreed upon. The first is that the Solar System and 
the low density clouds around it are embedded within a large, 50 — 200 pc 
radius, cavity that has a very low density of both neutral and ionized gas. Even 
allowing for large thermal pressure variations in the ISM, we expect that this 
volume is filled with a very low density ionized plasma, which provides some 
degree of pressure support for the cavity. The view that held sway for many 
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years was that the Local Bubble was filled with hot, relatively high pressure 
gas with a temperature of ^ 10^ K and pressure of P/k-Q ~ 10"^ cm~^ K. 
These conclusions were drawn from observations of the diffuse soft X-ray 
background (SXRB), which was mapped over the entire sky at low spectral 
and spatial resolution but high sensitivity at energies from ~ 70 eV to several 
keV by the Wisconsin group using sounding rockets (e.g. McCammon et al. 
1983). 

The ROSAT all-sky survey map, which has higher spatial resolution and 
sensitivity than the Wisconsin group maps, has made it clear that a substantial 
amount of the SXRB at high latitude comes from beyond the boundaries of the 
Local Bubble. X-ray shadows made by clouds outside the Local Bubble have 
provided information about the fraction of the emission coming from within 
the cavity. 

Another complication regarding the source of the soft X-ray emission is that 
it may be contaminated by X-ray emission from gas much closer to Eaith. As 
discussed by Cox (1998), and further explored by Cravens (2000), Cravens 
et al. (2001), and Pepino et al. (2004), a potentially substantial contribu- 
tion to the SXRB may be coming from the heliosphere and geocorona. The 
emission mechanism involves charge transfer reactions of highly charged ions 
such 0+^ in the solar wind with neutral H or He either in the interstellar wind 
or geocorona. After the charge transfer, the ion is in a highly excited state 
and radiatively decays to ground by emitting Lyman or K-shell photons. The 
heliospheric emission is expected to peak in emissivity at about 1 — 10 AU 
because of the combination of the spatial dependence of the ionization of in- 
flowing neutral atoms and the density of the solar wind. Current estimates are 
still quite uncertain, but put the contribution of the charge transfer photons to 
the observed SXRB at about 10 - 50% (Wai-gelin et al., 2004, Cravens et al., 
2001). 

It is important to note that even if the heliospheric and geocoronal soft X-ray 
emission is at the upper end of current estimates, they are unimportant for LIC 
ionization because the flux in the LIC, which is much further away ( ^ 100 x), 
is insignificant. Moreover, the emission observed by ROSAT and the Wisconsin 
sounding rockets ( ^ 70 eV) is substantially harder than the emission that is 
directly responsible for LIC ionization (~ 13.6 — 54.4 eV). Thus the observed 
SXRB is relevant to the photoionization of the cloud only insofar as it provides 
information on the hot gas that is emitting EUV photons that ionize the LIC. 

The available information on the temperature of the hot gas in the Local 
Bubble, which we still believe to be responsible for most of the observed 
SXRB, particularly near the Galactic plane, is mostly at low spectral resolution. 
The all sky maps of ROSAT in the energy band near 1/4 keV have been used 
to infer an effective plasma temperature near 10^ K (e.g. Kuntz and Snowden, 
2000) under the assumption of a plasma in coUisional ionization equilibrium 
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(CIE). The two higher spectral resolution data sets that are available (DXS, 
Sanders et al., 2001 and XQC, McCammon et al., 2002) are for very limited 
portions of the sky and at very low spatial resolution, and are not consistent 
with this model. In fact these spectra do not fit any currently known model 
for the temperature and ionization state of the plasma. One thing that does 
seem clear at this time is that iron must be significantly depleted in the plasma 
because the bright Fe line complex near 70 eV is observed to be much fainter 
than predicted. Perhaps this should not be surprising, since in the LIC ^^95% 
of the Fe is missing from the gas and presumably depleted onto dust grains, 
while in cold clouds ~99% of the Fe is depleted. 

Disregarding the problems with models of the emission spectrum, under the 
assumption that there is a CIE hot plasma filling the LB, the implied density is 
roughly 5 x 10"^ cm~^ and pressure of more than P/kB = 10"^ cm~^K. This 
is a high thermal pressure for the ISM, though certainly not unheard of, and 
may be typical of hot gas in the ISM. If the LIC is only supported by thermal 
pressure, however, there would seem to be a substantial mismatch between its 
pressure, P/kB ~ 2300 cm-^K, and that of the LB (§6.3.4). This mismatch 
may be fixed by magnetic pressure, though the size of the heliopause appears 
to require a field somewhat below that necessary to make up the difference 
(Ratkiewicz et al., 1998). That limit can also be avoided for a magnetic field 
that is nearly aligned with the direction of the interstellar wind (Florinski et al., 
2004), since in this case the heliosphere confinement is much less affected 
by the field. If a substantial portion of the SXRB originally attributed to the 
hot gas is from charge transfer emission, then models that properly take that 
into account will necessarily imply lower pressures for the LB, which helps 
to resolve this discrepancy. Note, though, that the emissivity goes as P^, so a 
50% reduction in emission results in only a factor of 1/ \/2 reduction in P. 

6.3 Neighborhood ISM: Cluster of Local Interstellar 
Clouds 

The ISM surrounding the Sun is flowing through the LSR with a bulk ve- 
locity of —19.4 lb 4.6 km s~^ and an upstream direction , —5°, 
towards the Scorpius-Centaurus Association (SCA). This warm, low density 
gas consists of cloudlets defined by velocity, which we denote the cluster of 
local interstellar cloudlets, CLIC. The upwind direction of the bulk velocity of 
the CLIC in the LSR suggests an origin related to stellar activity in the SCA 
(Frisch, 1981, Frisch and York, 1986, Frisch, 1995, also see §6.6.4). Low col- 
umn densities are found for ISM within 30 pc, Af(H°)<10^'^ cm^^. The CLIC 
temperature, iV(H°), and kinematics resemble global warm neutral medium 
(WNM) detected by radio H° 21 cm data (§6.6). 
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If we restrict the star sample to objects within ~10.5 pc, the mean A^(H°) 
for cue cloudlets is (iV(H°)) ~ 6.4 x 10^^ cm^^ based on H° or D° UV 
data (§6.1.1, and Table 6.2). For n(H°) similar to the LIC density at the solar 
location, n(H°)=0.2 cm^^, the mean cloud length is 1.0 pc, and the mean time 
for the Sun to cross them is ~68,000 years (Table 6.2). The LSR solar velocity 
is 13 — 19 km s~^ (§6.1.2), and the solar galactic environment over timescales 
of 3 Myrs will be regulated by ISM now within ^^60 pc, provided there are no 
undiscovered high velocity {V >10 km s~^) clouds with suitable trajectories. 



Table 6.2. The Characteristics of ISM within 10 pc of the Sun 



LSR Upwind direction 


£~331.0°,6~-5.r 


LSR Velocity(^) 


l/=-19.4±4.6 km s"i 


No. stars sampling ISM within 10 pc^'^) 


20 


SightUne averaged H° space density 


(ra(H°))= 0.07 cm-3 


No. velocity components''^) 


30 


Component averaged, A^(H°) 


(7V(H°))= 6.4 X IQi'^ cm-2 


Component averaged cloud thickness, L, 


(L) = 1.0 pc 


for n{W)=0.1 cm-3 




Component averaged cloud crossing 


(Tl) = 68,000 years 


time, Tl 




Galactic center hemisphere filling factor, / 


0.40 


Anti-center hemisphere filling factor, / 


0.26 


Solar entry into the CLIC 


(44,000-140,000)/nH/,o.2 yra 


Solar entry into the LIC 


<40,000/nH/,o.2 yra 


Solar exit from LIC 


next ~3700/n///,o.2 years 



Notes: The unit "years ago" is represented as "yra". nHi,o.2 is the average H° 
density in units of 0.2 cm^'^. (a) From Table 6.1. (b) The HD numbers of these 
stars within 10 pc are: 10700, 17925, 20630, 22049, 23249, 26965, 39587, 
48915, 61421, 62509, 48915B, 115617, 128620, 128621, 131156, 155886, 
165341, 187642, 197481, 201091, 209100. Sources of data for these stars are 
given in §6.1.1. (c) Based on UV data with resolution ~3 km s^^, generally. 

The first evidence of interstellar gas within 15 pc was an anomalously strong 
Ca+ line observed towai^ds Rasalhague (a Oph), formed in what is now known 
as the G-cloud that is widespread in the galactic-center hemisphere (Adams, 
1949, Munch and Unsold, 1962). The first spectral data of Lya emission from 
interstellar gas inside of the heliosphere, obtained by Copernicus, showed the 
similarity of ISM velocities inside of the heliosphere and towards nearby stars 
(Adams and Frisch, 1977, McClintock et al., 1978). 

The basic properties of ISM forming the galactic environment of the Sun 
were discovered with Copernicus, including the widespread presence of par- 
tially ionized gas (§6.6.1), the asymmetrical distribution of local ISM showing 
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higher column densities towards the galactic center hemisphere (Bruhweiler 
and Kondo, 1982, Frisch and York, 1983), and the discrepancy between the 
velocity of ISM inside of the solar system and towards the nearest star a Cen 
(Landsman et al., 1984, Adams and Frisch, 1977). The shocked history of the 
nearest ISM was revealed by enhanced abundances of Fe and other refractory 
elements, which indicated the destruction of interstellar dust grains by inter- 
stellar shocks (Snow and Meyers, 1979, Frisch, 1979, Frisch, 1981, Crutcher, 
1982, York, 1983). The first data showing the shift between the velocities of 
H° inside of the heliosphere and towards nearest star a Cen, now interpreted 
as partially due to the hydrogen wall, were obtained by Copernicus and lUE 
(Landsman et al., 1984), although the hydrogen wall contribution was not rec- 
ognized as such until the models of H° deceleration in the heliosheath were 
constructed (Gayley et al., 1997, Linsky and Wood, 1996). 

Optical or UV data are now available for ~100 stars sampling nearby ISM 
(see Frisch et al., 2002, Redfield and Linsky, 2004b, Wood et al., 2005, and ref- 
erences in these papers). Very high-resolution optical data, ~0.3-0.5 km s~^, 
are available for ^^40 nearby stars, and high resolution UV data, ~3 km s~^ for 
an additional ^■65 stars. Component blending between local and distant ISM 
usually prevent the use of data from distant stars. The H° La line is always sat- 
urated, even for low column density sightlines, so we use D° as a proxy for H°, 
with a ratio D°/H°=1.5 xlQ-^ that is valid for local ISM (Vidal-Madjar and 
Ferlet, 2002, Linsky, 2003, Lehner et al., 2003). These data, combined with ra- 
diative transfer models of ionization gradients, give the velocity, composition, 
temperature, and morphology of the CLIC. 

6.3.1 Warm Partially Ionized Medium, WPIM 

Charged and neutral particles have different interactions with the helio- 
sphere, and the ionization gradient of the cloud affects the heUosphere as it 
traverses a tenuous cloud. Over 30 years ago Copernicus discovered that N"'" 
is widespread towards nearby stars, N"'^/N°~l, indicating that EUV photons 
capable of ionizing hydrogen, nitrogen, and oxygen penetrate to the interi- 
ors of tenuous clouds and showing that the nearest ISM is partially ionized, 
X(H)>0.1 (Rogerson et al., 1973). This discovery contradicted the classic view 
of fully neutral or fully ionized ISM derived from observations of denser ISM, 
log A^(H°) (cm~^) > 20, where electrons originate from low FIP, < 13.6 eV, 
abundant elements such as C. The ionizations of H, O, and N, with FIP = 13.6, 
13.6, and 14.5 eV respectively, are coupled by charge exchange. The excep- 
tion is near a cloud edge where high EUV flux causes N°/H° to dip because 
photoionization dominates charge exchange ionization for N° (§6.4, Fig. 6.8). 
Partially ionized gas in the CLIC is established by N fractional ionizations of 
X(N)~0.27-0.67, derived from FUV and UV data on N° and N+ towards stars 
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with logiV(H°)(cm-2) < 19, such as HZ 43, HD 149499B, WD 0549+158, 
WD 221 1-495, and r/ UMa (Lehner et al., 2003, Frisch et al, 2005). 

Among these stars the most highly ionized sightlines are HD 149499B and 
WD221 1-495. Both of these hot stars have detected O^^ absorption, which 
appears to be interstellar, although a stellai" origin can not be ruled out entirely 
(Oegerle et al., 2005). These stars are in the upwind hemisphere of the CLIC 
(§6.3.2), with HD 149499B at £=330°, h=-T , d=37 pc, and WD2211-495 at 
£=346°, 6=-53°, d=53 pc. In contrast, ISM towards WD 1615-154 is primar- 
ily neutral with x(N)~0.05 (£=359°, 6=-i-24°, d=55 pc). Several stars sample 
mainly the LIC, such as WD0549-I-158 near the downwind direction (£=192°, 
6=-5.3°, d=49 pc), where iV(H°)~ 5 x 10^^ cm"^ is comparable to the ex- 
pected LIC column density. 

The H and He ionizations are found from EUV spectra of white dwarf stars, 
and comparisons of fluxes at the He° and He+ ionization edges at 504 A and 
229 A with atmosphere models. Well observed stars such as HZ 43, GD 153, 
and WD 0549-^158 show iV(H°)/iV(He°)~ 9.8 - 15.8, instead of the ratio of 
10/1 expected for neutral gas with a H/He = 10 cosmic abundance. When 
lower quality data are included, the variation is larger, A^(H°)/iV(He°)=9-40 
(Dupuis et al., 1995, Kimble et al., 1993, Frisch, 1995, Vallerga, 1996, Wolff 
et al., 1999). Ratios of H°/He° > 10 indicate that He is more highly ionized 
than H, and the hardness of the radiation field implied by this discovery is 
discussed in §6.4.1. 

The ratios iV(Mg+)/iV(Mg°) and iV(C+)/iV(C+*) serve as ionization di- 
agnostics and give n(e) values that are independent of abundance uncertain- 
ties. The Mg° abundance is enhanced by dielectronic recombination in warm 
ISM, T > 6000 K, and Mg'^ is the dominant ionization state. The elec- 
tron densities are given by n{e)=CMg(T, A) iV(Mg°)/Mg+ and n(e)=Cc(T) 
A^(C+*)/A^(C+). The quantity CMg(r,A) is the ratio of the temperature de- 
pendent photoionization and recombination rates, and depends on both cloud 
temperature (T) and the FUV radiation field (A). Cc(T) is the temperature 
sensitive ratio of the fine-structure coUisional de-excitation to excitation 
rates and is independent of the radiation field. For more details, see York and 
Kinahan (1979). The observed ratios in the CLIC of iV(C+)/iV(C+*)=50-200 
and A^(Mg'^)/A^(Mg°)=200^50 are consistent with the predictions of equi- 
librium radiative transfer models of low column density gas such as the LIC 
(section §6.4). These equilibrium models predict n(e)=0.06-0.12 cm"'^, with 
the best value for the CLIC being '^O.l cm"'^. Some caution is required when 
using observed C"*" column densities, however, since the 1335 A line is 
generally saturated, which may bias A^(C"*")/A^(C"*"*) towards smaller values. 

Neutral Ar, FIP=15.8 eV, is a valuable ionization diagnostic since it is ob- 
served inside the solar system, in the form of the anomalous cosmic ray com- 
ponent seeded by interstellar- Ar°, and also in FUV observations of nearby 
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Stars. The ratio log Ar°/H° in the ISM, ~ —5.8 ± 0.3 dex, is below values 
seen towards B-stars by factors of two or more (Jenkins et al., 2000, Lehner 
et al., 2003). The recombination coefficients of Ar° and H° are similar, but the 
Ar° photoionization cross section exceeds that of H° by ~10. Argon deple- 
tion is expected to be minimal, so Ar/H towards nearby stars is interpreted as 
conversion to unobserved Ar"*". The RT models predict Ar+/Ar° = 1.8 — 2.9, 
Ar°/Ar = 0.16 - 0.30, and logAr°/H° = -5.98 to - 6.17 for the LIC 
at the solar location, compared to observed ACR values of logAr°/H° = 
—5.97 ± 0.16 after filtration effects in heliosheath regions are included (see 
FS for more discussion and references). 

A classic ionization diagnostic uses the strong optical lines of the trace ion- 
ization species Na° and Ca^and the assumption that ionization and recombi- 
nation rates balance: 

n(Na°)r(Na°) = n(Na+)Q(Na+)n(e) (3.1) 

Implicit is the assumption that A^(Na°)=L n(Na°) cm~^, or analogously that 
n(Na°)/n(Na+) is constant, where L is the cloud length. r(Na°) is the total 
ionization rate for Na°— >Na+, and a(Na°) is the total recombination rate for 
Na+^Na° (e.ge.g. Spitzer, 1978, Pottasch, 1972). In principle Ca+ andCa++ 
can be substituted for Na°, and Na+, respectively. Na"*" and Ca"'"^" are unob- 
servable. Since Na° (and Ca+) are trace ionization states in warm gas, some 
assumption is required for the element abundances, (and [^]): 

n(Na°)r(Na°) = n(H)[^]a(Na+)n(e) (3.2) 

H 

The electron densities derived from Na° and Ca"*" data may be unreliable for 
low column density warm clouds because of large ionization corrections and 
highly variable abundances, of ~40 and ~1.6, respectively, for Ca and Na 
(Welty et al., 1999, §6.1.1). Calcium ionization levels are temperature sen- 
sitive, and Ca++/Ca+>1 for T > 4,000 K and n(e)<0.13 cm-^. Electron 
densities of n(e)=0.04-0.19 cm~^ are found for the CLIC within 30 pc from 
Mg+/Mg° and Na° (Frisch et al., 1990, Lallement and Ferlet, 1997). 

The ratio Fe+/D° presents clear- evidence for variations in the ISM over spa- 
tial scales of 1-3 pc, however, it is unclear whether the variation arises from 
an ionization gradient or from a gradient in the Fe abundance close to the Sun. 
Ionization and depletion both affect Fe^/D°, since Fe'^ is the dominant ion- 
ization state of Fe for both neutral and warm ionized gas (SF02) and D° traces 
only neutral gas. Fig. 6.3 shows Fe+/D° for CLIC velocity components, plot- 
ted as a function of the angle, 6, between the star and the LSR upwind direction 
of the CLIC (Table 6.1). The ratio iV(Fe+)/iV(D°) for cloudlets within ~4 pc 
differs by ~54% over spatial scales of 4 pc between upwind stai^s such as a Cen 
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AB, ~16°, and downwind stars such as Sirius, d ~ 103°, although uncer- 
tainties overlap. Similar Fe abundance variations are found in thin (~ 0.1 — 0.9 
pc) disk clouds towards Orion (Welty et al., 1999). The Fe abundance is highly 
variable in the ISM, and varies between cold dense and warm tenuous disk 
clouds by factors of ~4-6 (Savage and Sembach, 1996, Welty et al., 1999). 
The local Fe+/D° trend may be from a combination of ionization and abun- 
dance variations. The first explanation for the physical properties of local ISM 
attributed the CLIC to a superbubble shell around the SCA that is now ap- 
proaching the Sun from the upwind direction (Frisch, 1981, §6.6.4). Dust grain 
destruction in the shocked superbubble shell would have restored elements in 
silicates, which can be destroyed by shocks with speeds > 50 km s~^, back 
to the gas phase (Slavin et al., 2004, Jones et al., 1994). This effect is echoed 
in extended distances to the upwind edge of the CLIC obtained from Ca+ data 
(Figure 6.5). 

A diffuse H II region appears to be close to the Sun in the LSR upwind direc- 
tion of the CLIC, based on ISM towards A Sco, and the high value x(N)~0.66 
towards HD 149499B (WD 1634-593, Lehner et al., 2003). About 98% of the 
neutral gas towards A Sco belongs to the CLIC, and shows a HC velocity of 
— 26.6 km s~^. A diffuse H II region, with a diameter of ^^30° and n(e)~0.1- 
0.3 cm~^, is also present at ~ —17.6 km s~^ (using optical data to correct 
Copernicus velocities to the HC scale, York, 1983). This diffuse HII region 
may extend in front of both the HD 149499B and A Sco sightiines. The HD 
149499B sightline samples the ISM approaching the solar system from the 
LSR CLIC upwind direction. 

Highly ionized gas is observed in the downwind CLIC. Interstellar Si"*"^ is 
seen towards e CMa but not towards a CMa (Sirius), ~12° away (Gry and 
Jenkins, 2001, Hebrard et al., 1999). The limit on Si''"^ towards Sirius is a fac- 
tor of 15 below the e CMa LIC component, with iV(Si+2) = 2.3 ± 0.2 x 10^^ 
cm~^. This ion is particularly interesting because it is predicted to have quite a 
low column density both in the warm LIC gas, because of its high charge trans- 
fer rate with H°, and in an evaporative boundary where it becomes quickly 
ionized in the interface (see SF02). Gry and Jenkins speculate that Si"*"^ is 
formed in an outer layer of the LIC behind Sirius (2.7 pc). However this sce- 
nario requires that the second local cloud observed towards both stars, which is 
blue-shifted by ~ —7 km s~^ from the LIC (hence the name Blue Cloud, BC), 
is a clump embedded in an extended LIC. Since Si+^ is also detected for the 
BC towards e CMa but not Sirius (but with lai^ge uncertainties), the origin of 
the Si+^ towards e CMa may instead require nonuniform interface layers such 
as turbulent mixing layers (see below). 
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Figure 6.3. The ratio A^(Fe^)/A'^(D°) varies between ttie upwind and downwind direction 
of the CLIC, and over spatial scales of ~ 3 pc. Ai'(Fe+)/A'^(D°) is plotted against the an- 
gle (in degrees) between the background star and the LSR upwind direction of the CLIC 
{£ = 331.4°, & = —4.9°, Table 6.1). Clouds with observed velocities within ±1 km s~^ of 
the projected LIC velocity are plotted as filled circles, while open circles represent CLIC com- 
ponents at other velocities. The open squares represent lower resolution data towards a Oph, A 
Sco and HD 149499B. The ratio Af(Fe+)/iV(D°) for the LIC differs between the direction of 
a Cen AB, at a distance 1.3 pc and angle of 16°, and Sirius, at a distance of 2.7 pc and angle 
of 103°, although uncertainties overlap. Velocities are drawn from the data compilations in the 
RL and Frisch et al. (2002) surveys. 



6.3.2 Dynamical Characteristics of Nearby ISM 

The heliosphere radius in the upwind direction is approximately propor- 
tional to the relative velocity between the Sun and ISM, so that the heliosphere 
is modified by variations in the ISM velocity. Absorption line data give the 
centroid of the radial component of a cloud velocity, integrated over a cloud 
length, but UV data do not resolve all of the velocity structure (§6.1.3). The 
coherent motion of nearby ISM is found from these velocity centroids. With 
the exception of the LIC, which is detected inside of the heliosphere, only ra- 
dial velocities can be measured. Thus the 3D motions of other cloudlets must 
be infeiTcd from observations towards several stars. 

It has been known for some time that the Sun is immersed in ISM flowing 
away from the Scorpius-Centaurus Association and towards the Sun (Frisch, 
1981, Crutcher, 1982). The CLIC kinematical data can be interpreted as a co- 
herent flow (Frisch and York, 1986, Vallerga et al., 1993), or individual clouds 
can be identified with similar velocities towards adjacent stars (Lallement et al., 
1986, Frisch et al., 2002). Clouds that have been identified are listed in Table 
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6.3. Among these clouds is the G-cloud, in the galactic center hemisphere, 
which was first identified in optical data long ago (Adams, 1949, Munch and 
Unsold, 1962). The Sun is located at the leading edge of the stream of CLIC 
gas (Frisch, 1995). 

Kinematics of nearby ISM show three characteristics. The first is that the 
ISM within ~30 pc flows past the Sun. The bulk flow velocity, Vclic> can 
be derived from the Doppler-shifted radial velocities of ~100 interstellar ab- 
sorption line components towards ~70 nearby stars (Frisch et al., 2002). The 
resulting Vclic corresponds to a heliocentric velocity -28.1 ±4.6 km s~^ and 
upwind direction (^,6)=(12.4°,1 1.6°) (Frisch et al., 2002). An additional un- 
certainty of ~ 1 km s~^ may be introduced by a biased sample, since there 
are more high-resolution optical data for stars in the upwind than in the down- 
wind directions. The CLIC motion in the LSR is then found by subtracting 
the solar apex motion from the CLIC heliocentric vector. Table 6. 1 gives the 
LSR CLIC velocity for two choices of the solar apex motion. Fig. 6.4 (top) 
displays the velocities of individual absorption components in the rest frame 
of the CLIC versus the rest frame of the LSR. Obviously the CLIC velocity is 
more representative of the nearby gas. 

The second property is that distinct clouds in the flow contribute to the ±4.6 
km s~^ dispersion of Vclic (Tables 6.3 and 6.1). Among the clouds are the 
Apex Cloud, within 5 pc and centered towards the direction of solar apex mo- 
tion, and the G-cloud, which is seen towards many stars in the upwind hemi- 
sphere. The Apex and G-cloud data indicate that temperature varies by a factor 
of 2-6 inside of the clouds (Table 6.3, §6.3.4, RL). Evidently, for constant pres- 
sure the clouds are either clumpy, or there is a dynamically significant magnetic 
field. For densities similar to the LIC, n(H°)~0.2 cm~^, the clouds fill only 
~ 33% of nearby space and the flow of warm gas is fragmented, rather than 
a streaming turbulent homogeneous medium (§6.3.3). The interpretation of 
CLIC kinematics as a fragmented flow is supported by the fact that the upwind 
directions for the CLIC, LIC, Apex Cloud, and GC are all within 20° of each 
other. 

A third characteristic is that the flow appears to be decelerating. This is 
shown in Fig. 6.4, bottom. Stars in the upwind and downwind directions are 
those with projected CLIC velocities of < km s^^and > km s~^, respec- 
tively. The velocities of clouds closest to the upwind direction are approaching 
the Sun compared to Vclic> while those closest to the downwind direction lag 
Vclic > as would be expected from a decelerating flow. This deceleration is 
seen even for the nearest stars, d < 6 pc, indicating the pileup of ISM is close 
to the Sun, which is consistent with the fact the Sun is in the leading edge of 
the flow. 

Note that because the 3D LIC velocity vector is known, a 3D correction 
for the solar apex motion is made, and the "true" downwind direction in the 
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Figure 6.4. Top: The star longitude is plotted against the velocities, dV, of CLIC absorption 
components after transforming the absorption components into the Vclic frame (filled circles), 
and into the LSR (open circles. Table 6.1). Clearly Vclic provides a better fit to the CLIC data 
than the LSR velocity frame. LSRstd is used here, and only components within 50 pc of the 
Sun are plotted. Bottom: The observed heliocentric velocities, V, of absorption components 
tracing the CLIC are plotted against the expected CLIC velocity projected towards each star 
(Table 6.1). The components near the upwind (approaching gas, V <Q km s~^) and downwind 
(receeding gas, V > Q km s^^) directions, which are expected to have the maximum values of 
projected | V|, show indications of a decelerating flow. 



LSR differs from the direction measured with reference to the solar system 
barycenter. 

6.3.3 The Distribution of Nearby ISM 

The distribution of ISM within 5-40 pc is dominated by the CLIC, while 
over larger scales it is dominated by the Local Bubble (§6.2.1). The >15° 
diameter cloud observed towards the Hyades stars, ^■40-45 pc away, is not 
included in the CLIC discussion here because A^(H°) data aie unavailable, 
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Table 6.3. Nearby Interstellar Clouds 



Cloud 


Velocity 


Location 


Temp. 




HC 


LSR 


D,ib 






(kms-i) 


(kms-i) 


(pc,° °) 


(10^ K) 


Lie 


-26.3±0.3 


-20.7 


0., 50° ^250° (all b) 


6.4±0.3 


G 


-29.3±4.0 


-21.7 


1.4, 0°±90° 0°±60° 


2.7-8.7 


Blue 


lOibl 




3, 233° ±7° 10°±2° 


2.0-5.0 


Apex 


-35.1±0.6 


-24.5 


5, 38°±10° 9°±14° 


1.7-13.0 


Peg/Aqr 


-4. 5 ±0.5 




30, 75°±13° ^4° ±5° 





Notes: The LSR velocities are based on the Standard solar apex motion. The 
cloud distances are upper limits. The Apex cloud is also known as 
"Panoramix" or the "Aql-Oph cloud". References: Table 6.1, Witte (2004), 
FGW, Frisch (2003), Lallement et al. (1986), Lallement and Bertin (1992),GJ, 
RL. 



and because the nonlocal ISM may be inside the cluster (Redfield and Linsky, 
2001). 

The solar apex motion is compared to the distribution of CLIC gas in Fig. 
6.5, where the distance to the cloud edge for the CLIC is shown projected onto 
the galactic plane, and in a vertical plane perpendicular to the galactic plane 
and aligned along the solar apex motion. The angular width of the vertical dis- 
play, which extends from fc40°±25° to ^=220° ±25°, includes the directions 
of both the Hipparcos and Standard solar apex motions (Table 6. 1). The down- 
wind direction of Vclic and direction of solar motion (Table 6. 1) are shown by 
the arrows. The distance of the CLIC edge in the direction of a star is given by 
A^(H°)/(n(H°)), where A^(H°) is the total interstellar column density towards 
the star, and a uniform density similar to LIC values is assumed, (n(H°))=0.2 
cm"'^. This distribution is based on D°, H°, (dots) and Ca+ (crosses) data. 
The high resolution optical Ca"*" data provide excellent velocity resolution, but 
conversion to A^(H°) is uncertain because of possible variable abundances. The 
value A^(Ca+)/A^(H°)=10~^ cm~^ is used. The solar apex motion is shown by 
the arrows pointing right, while the arrows pointing left show the two LSR 
velocity vectors of the CLIC bulk flow, for the Standard (solid) and Hipparcos 
(dotted) solar apex motions, respectively. 

The CLIC shape is based on data in Hebrard et al. (1999), RL, Dunkin and 
Crawford (1999), Crawford et al. (1998), FGW, and Frisch and Welty (2005). 
The iV(H°) column densities are either estimated from iV(D°) (D°/H°= 1.5 x 
10~^), or based on the saturated H° La line, or estimated from A^(Ca"'") using 
iV(Ca+)/iV(H°)= lO"*^. The iV(Ca+)/iV(H°) conversion factor is based on 
CLIC data towards nearby stars such as a Aql, rj UMa, a CMa, and is uncertain 
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Figure 6.5. See figures fig5a.pdf fig5b.pdf. The distance to the CLIC edge for n(H°)'~0.2 
cm~^ and a continuously distributed ISM. The dots give distances derived from A'^(D°) and 
iV(H°) data, and the x's show Ca+ distances. The extended HII region found towards A Sco 
(§6.3. 1, York, 1983) is shown as the origin of the excess cloud length towards the galactic center 
indicated by the Ca^ data. It is interpreted to indicate fully ionized gas near the Sun because the 
diffuse H II region seen towards A Sco is at CLIC velocities. The arrows give the motions based 
on the standard (solid) and Hipparcos (dashed) solar apex motions for both the Sun and CLIC. 
Top: The edges of the CLIC are shown projected onto the galactic plane. Bottom: The CLIC 
distribution is shown for a meridian cut 50° wide in longitude, and extending between /=40° 
and 220° (the plane of the solar apex motion). The CLIC LSR motion is nearly perpendicular 
to the solar apex motion. 



because Ca depletion varies strongly. Most Ca is Ca"*""*" for warm gas, if T > 
4, 000 K and n(e)<0.13 cm~'^. Radiative transfer models of the LIC predict 
that Ca"'""'~/Ca~'"~40-50, so small temperature uncertainties may produce large 
variations in Ca+ZH". 

With the possible exceptions of a Oph, interstellar column densities for stars 
within 30 pc are less than log N{W){cvar'^) ~ 18.5 (Frisch et al., 1987, Wood 
et al., 2000a). Neutral CLIC gas does not fill the sightline to any nearby star if 
the ISM density is n(H°)~0.2 cm~'^. The CLIC extends farther in the galactic 
center hemisphere than the anti-center hemisphere, as traced by D° and Ca+ 
data. The most puzzling sightline is towards a Oph (14 pc), where Ca+ is 
anomalously strong and the 21 -cm H° emission feature at the same velocity 
suggests log A^(Htot)(cm~^) ~ 19.5 dex (see 6.5.3). High column densities 
(logiV(H°)(cm-2) ~ 18.75) are also seen towards HD 149499B, 37 pc away 
in the LSR upwind CLIC direction, and towards LQ Hya, where strong ISM 
attenuation of a stellar H° Lq emission feature causes the poor definition of 
the line core and wings. 

The percentage of a sightline filled with ISM offers insight into the ISM 
character, and is given by the filling factor, /. Restricting the discussion of the 
ISM filling factor to the nearest 10 pc, we find that ~67% of space may be 
devoid of H°. If n(H°)~0.2 cm'^, then /~0.33 for ISM within 10.5 pc. For 
galactic center and anti-center hemispheres, respectively, /~0.40 and /~0.26. 
Mean cloud lengths are similar for both hemispheres. The highest values are 
/~0.57, towards a Aql (5 pc) and 61 CygA (3.5 pc), and the lowest values 
are /~0.11 towards Ori (8-7 pc), which is ~14° from the downwind di- 
rection. The sightline towards Sirius (2.7 pc, 43° from the downwind direc- 
tion) has /~0.26. These filling factors indicate that the neutral ISM does not 
fill the sightline towards any of the nearest stars, including towards a Cen 
where /~0.5, unless instead the true value for (n(H°)) is much smaller than 
0.2 cm-3. 
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6.3.4 Cloud Temperature, Turbulence, and Implications 
for Magnetic Pressure 

The basic thermal properties of the CLIC are presented in the Redfield and 
Linsky (2004) survey of ~ 50 cloudlets towards 29 stars with distances 1 — 95 
pc. Cloudlet temperatures (T) and turbulence (0 values ai^e found to be in the 
range T = 1, 000 - 13, 000 K and ^=0-5.5 km s'^ for clouds within ~100 pc 
of the Sun. The mean temperature is 6680±1490 K, and the mean turbulent ve- 
locity is 2.24ibl.03 km s~^. From these values, RL estimate the mean thermal 
(PTh/k = nT) and mean turbulent (P^/k^O-Sp^^/j,) pressures of 2,280±520 K 
cm~^ and 89±82 K cm~^, respectively, by assuming n(H°)= 0.1 cm~^ and 
n(e)=0.11 cm"'^. The thermal pressure calculation includes contributions by 
H°, He°, electrons, protons, and assumes that He is entirely neutral. The pres- 
sure will be underestimated by ~15% if He is ^^50% ionized as indicated by 
radiative transfer models (§6.4). For comparison, if these clouds have ioniza- 
tion similar to the LIC, then PTh~2300 cm"^ K. 

These cloud temperatures are determined from the mass dependence of line 
broadening using the Doppler parameter, 6d> so spectral data on atoms or ions 
with a large spread in atomic masses are needed. In practice, observations 
of the D° La line are required for an effective temperature determination that 
distinguishes between thermal and nonthermal broadening. 

When the star sample is restricted to objects within 10.5 pc, the cloud tem- 
perature is found to be anticoiTclated with turbulence, and to be correlated 
with A^(D°) (Fig. 6.6). From a larger sample of components, RL have con- 
cluded that the T-^ anti-con^elation is significant. However the likelihood that 
unresolved velocity structure is present in these UV data allows for the T-^ 
anti-correlation to contain some contribution from systematic eiTors. High res- 
olution optical data show that velocity crowding for interstellar- Maxwellian 
components persists down to component separations below 1 km s~^ (§6.1.3), 
so that the weak positive correlation between T and A^(D°), and negative cor- 
relation between T and ^ may result from unresolved component structure. 

There are no direct measures of the magnetic field strength in the LIC, but 
the field strength is presumed to be non-zero based on observations of polar- 
ized starlight for nearby stars, which may originate from magnetically aligned 
grains trapped in interstellar magnetic field lines draped over the heliosphere 
(Frisch, 2005). The thermal properties of the CLIC have implications for pres- 
sure equilibrium and magnetic field strength. The magnetic field strength and 
density fluctuations can be constrained using equipartition of energy argu- 
ments. If clouds in the CLIC are in thermal pressure equilibrium with each 
other, P'Yh/k = nT, where n is the total number of neutral and charged par- 
ticles in the gas, and if magnetic field i? = 0, then the temperature range 
of T ~ 10^ — 10^ K found by RL indicates that densities must vary by an 
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Figure 6.6. Temperature versus turbulence ^ (right) and temperature versus column density 
A'^(D°) (left ) for interstellar absorption components seen towards stars within 10 pc of the Sun 
(based on data in RL). 



order of magnitude. Since particle number densities vary by a factor of ~2 
as the cloud becomes completely ionized, most of the temperature variation 
must be balanced either by variations in the mass-density or in the magnetic 
field strength if the cloud is in equilibrium. If the CLIC has a uniform total 
density, and if thermal pressure variations are balanced by magnetic pressure, 
Pe/fc = B'^/Skir, then magnetic field strengths must vary by factors of ~3 in 
the CLIC. If the star set is restricted to objects within ~10 pc, a temperature 
and turbulence range of T=l,700-12,600 K and ^ = — 5.5 km s^^ are found, 
with mean values of 6, 740 ± 2800 K and ^ = 1.9 ± 1.0 km s'^. 

A rough estimate is obtained for the magnetic field strength in the LIC by 
assuming equipartition between thermal and magnetic energies, and using the 
results of the RT models that predict neutral and ion densities (see §6.4). The 
first generation of models gives a LIC thermal energy density of E^h/k ~ 3600 
cm-3 K for T = 6340 K, and including H°, H+, e", He°, and He+. Lower 
total densities and ionization in the second generation of models reduce the 
thermal energy density somewhat. Equipartition between thermal and mag- 
netic energy density gives £'Th=^B and E^/k = B'^ /Skir. These assumptions 
then give S~ 3.1 - 3.8 fiG for the LIC. 

The interstellar magnetic field strength in the more extended CLIC can be 
guessed using the RL value for the mean thermal pressure of 2280 cm"'^ K, and 
assuming that the mean magnetic and thermal pressures are equal. For this case 
Br^l.S fiG. If these clouds are, instead, in pressure equilibrium with the Local 
Bubble plasma, then (Pxh+^B/)^ = PhB/k ~ 5 x 10^ K cm"^ (§6.2.2), and 
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magnetic field strengths are ~3.1 fiG. In contrast, for Pb~^{ then 5~0.6 fiG. 
Based on equipartition of energy arguments, typical field strengths of B ~3 
/iG seem appropriate for the CLIC, with possible variations of a factor of 3. 

The Lie turbulence appears to be subsonic. Treating the LIC as a perfect 
gas, the isothermal sound speed is Vs'^ 0.09VT ~ 7.1 km s~^, and turbulent 
velocities are 0.5-2.7 km s"^ (Hebrard et al., 1999, Gry a nd Jenkins, 2001, 
RL). The Alfven velocity is given by Va'^ 2.2i?^/y^n(p), where Bj^i is the 
interstellar magnetic field in 10~^ G {fiG), Va is in km s~^, and the proton 
density n(p) is in cm"'^. For gas at the LIC temperature (6,300 K), the Alfven 
velocity exceeds the sound speed for > 1.3 /iG. The velocity of the Sun 
with respect to the LIC (26.3 km s~^) is both supersonic and super-Alfvenic 
for interstellar field strengths B < 3.7 /iG. 

6.4 Radiative Transfer Models of Local Partially Ionized 
Gas 

Radiative transfer (RT) effects dominate the ionization level of the tenu- 
ous ISM at the Sun. The solar environment is dominated by low opacity 
ISM, A^(H°) ^ 10^^-^ cm~^. In contrast to dense clouds where only pho- 
tons with A >912 A penetrate to the cloud interior, the low column density 
ISM near the Sun is partially opaque to H-ionizing photons and nearly trans- 
parent to He-ionizing photons. At 912 A the cloud optical depth r ~1 for 
logiV(H°)~ 17.2 cm~^, and at the He° ionization edge wavelength of 504 A, 
r ~1 for logA^(H°)~ 17.7 cm~^. The average H° column and mean space 
densities for stars within 10 pc of the Sun are {N{Yl°)) ~ lO'^^ cm~^ and 
{n{¥l°))r^ 0.07 cm~'^, so the heliosphere boundary conditions and the ratio 
n(H°)/n(H+) vary from radiative transfer effects alone as the Sun traverses 
the CLIC (Fig. 6.8). Warm, T > 5000 K, partially ionized gas is widespread 
near the Sun and is denoted WPIM (§6.3.1). Charged interstellar particles cou- 
ple to the interstellar magnetic field and are diverted around the heliopause, 
while coupling between interstellar neutrals and the solar wind becomes sig- 
nificant inside of the heliosphere itself. The density of charged particles in the 
ISM surrounding the Sun supplies an important constraint on the heliosphere, 
and this density varies with the radiation field at the solar location, which is 
now described. 

6.4.1 The Local Interstellar Radiation Field 

The interstellar radiation field is a key ingredient of cloud equilibrium and 
ionization at the solai^ position. This radiation field has four primary compo- 
nents: A. The FUV background, mainly from distant B stai^s, B. Stellar- EUV 
emission from sources including nearby white dwarfs and B stars (e CMa and 
P CMa); C. Diffuse EUV and soft X-ray emission from the Local Bubble hot 
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plasma, as we discussed in §6.2.2; and D. Additional diffuse EUV emission 
thought to originate in an interface between the warm LIC/CLIC gas and the 
Local Bubble hot plasma (§6.2.2). This last component is required because, al- 
though radiative transfer models show that the stellar EUV and Local Bubble 
emission account for most LIC ionization, it is not sufficient to account for the 
high He ionization inferred throughout the cloud from the EUVE white dwarf 
data (Cheng and Bruhweiler, 1990, Vallerga, 1996, Slavin, 1989). The spectra 
of these radiation sources are shown in Fig. 6.7. 

The interstellar radiation flux at the cloud surface must be inferred from 
data acquired at the solar location, together with models of radiative transfer 
effects. LISM column densities are so small that dust attenuation is minimal, 
e.g. for the LIC Ay < 10^^ mag, and fluxes longwards of A ~ 912 A are 
similar at the solar location and cloud surface (with the exception of Ly a 
absorption at 1215. VA). For wavelengths A < 912 A the situation is different, 
however, and the spectrum hardens as it traverses the cloud because of the high 
H°-ionizing efficiency of 800 — 912 A photons. Thus a self-consistent analysis 
is required to unravel cloud opacity effects, and extrapolate the EUV radiation 
field observed at the Sun to the cloud surface. The observational constraints 
on the 200 — 912 A radiation field are weak, partly due to uncertainties in 
A^(H°) towards e CMa and partly due to the difficulty in observing diffuse 
EUV emission, which allows some flexibility in introducing physical models 
of the cloud interface. 

A thin layer of intermediate temperature ISM is expected to exist in the 
boundary between the CLIC and the LB hot plasma. This layer will emit ra- 
diation with a spectrum and flux dependent on the underlying physical mech- 
anisms. Models for this interface emission indicate it radiates strongly in the 
energy band £^ = 20 — 35 eV, which is important for both He° and Ne° ion- 
ization. The exact physical processes at work in interface regions are unclear, 
but possibilities include thermal conduction, radiative cooling, and shear flow, 
which lead to evaporative interface boundaries (Cowie and McKee, 1977), 
cooling flows (Shapiro and Benjamin, 1991) or turbulent mixing layers (Slavin 
et al., 1993). All of these boundary types produce intermediate temperature 
gas (T ~ 10^'^ — 10^'^ K) that radiates in the EUV, although ionization levels, 
and thus the spectrum and intensity of the emission, depend on the detailed 
physics. Parameters that constrain interface properties include the strength and 
topology of the interstellar magnetic field. Bis, the hot gas temperature, and the 
relative dynamics of the hot and warm gas. The magnetic field affects the RT 
models because Bjs reduces the evaporative flow by inhibiting thermal con- 
duction in directions perpendicular to field lines and at the same time supports 
the cloud by magnetic pressure. Since the total (magnetic + thermal) pressure 
is roughly constant in an evaporative outflow, the magnitude of the magnetic 
field is an important factor in determining the pressure in any evaporative flow 
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that might be present. Fig. 6.7 shows examples of the EUV spectrum produced 
by an evaporative interface and a turbulent mixing layer. 




Figure 6.7. Left: Components of the interstellar radiation incident on the local interstellar 
cloud complex. Contributions from stars show the EUV flux from nearby white dwarf and B 
stars, after deabsorbing the corresponding Af(H°) to the cloud surface. The cloud boundary flux 
produced by emission from an evaporative interface between the local gas and the hot gas of 
the Local Bubble is shown (Slavin and Frisch, 2002). Right: Radiation field from a turbulent 
mixing layer interface compared to that of a conductive interface. Instabilities at the boundaries 
of tenuous clouds may lead to quite different radiation fields inside of the clouds, particularly 
near 584 A where He° is ionized. 



6.4.2 Radiative Transfer Models of the Local Cloud and 
other Tenuous ISM 

Our radiative transfer models are constrained by observations of local ISM 
towards nearby stars, and in situ data from LIC neutrals that have penetrated 
the heliosphere. The in situ data includes direct detection of He°, observations 
of solar La florescence from H° in the heliosphere, and observations of the 
pickup ion and anomalous cosmic ray populations that are seeded by interstel- 
lar neutrals (see Moebius et al.. Chapter 8). Generally elements with FIP< 13.6 
eV are fully ionized in tenuous clouds, while elements such as H, O, N, Ar, He, 
and Ne with 13.6 ^ PIP ^ 25 eV are partially ionized. Neutrals from these par- 
tially ionized species enter the heliosphere, where they seed the pickup ion and 
anomalous cosmic ray populations measured by instruments on various space- 
craft. Column densities towards nearby stars constrain sight-line integrated 
values, and permit the recovery of n(H'^)/n(H°) as a function of distance to 
the cloud surface (Fig. 6.9). While the LIC temperature is determined directly 
from Ulysses observations of He°, the densities n(H°) and n(H"*") at the solar 
location vary as the Sun moves through the LIC, and must be determined from 
radiative models. 



34 



Short-term Variations in the Galactic Environment of the Sun 



The detailed attention paid here to LIC radiative transfer models is moti- 
vated by the facts that ri(H°) and n(He°) at the solar location are important 
boundary conditions of the heliosphere, and that ionization gradients in the 
CLIC are factors in reconstructing the 3D cloud morphology from data. An 
extensive study of the LIC ionization has resulted in a series of ^^50 radia- 
tive transfer models appropriate for tenuous ISM such as the LIC and CLIC 
(Slavin and Frisch, 2002, Frisch and Slavin, 2003, Frisch and Slavin, 2005). 
These models use the CLOUDY code (Ferland et al., 1998), an interstellar ra- 
diation field based on known sources, and models for the interface expected 
between the LIC and Local Bubble plasma (also see §6.2.2). Boundary condi- 
tions for outer scales are provided by ISM data integrated over the LIC, which 
can be subject to ionization gradients. On inner scales the boundary condi- 
tions are given by in situ observations of ISM at the heliosphere. Both sets of 
constraints are important for evaluating the RT properties of the surrounding 
ISM where large H°/H+ gradients ai^e found. The data sets used here include 
n(He°), pickup ions, and anomalous cosmic rays inside of the heliosphere, and 
absorption line data of the LIC towards the downwind star e CMa (GJ, FS). 

These RT models will be updated in the future, as the density, composition, 
LIC magnetic field, radiation field, and interface regions between the LIC and 
Local Bubble plasma become better understood. The LIC magnetic field enters 
indirectly thi^ough the contribution of magnetic pressure to the cloud interface. 
The magnetic pressure in the cloud helps to determine the thermal pressure and 
thus cooling in the interface, so that further studies of i^is near the Sun will 
provide insight into the characteristics of the interface radiation. One of the 
significant results of our study is that the spread in predicted neutral and ion 
densities demonsti'ates that low column density ISM can be in equilibrium at 
a range of ionization levels, and that the ionization itself is highly sensitive to 
the radiation field, interface characteristics, and other cloud properties. 

6.4.3 Results of Radiative Transfer Models 

Two generations of radiative transfer (RT) models have been developed, 
with the focus on matching data for the LIC inside of the heliosphere, and 
matching data on nearby ISM in the downwind direction where the brightest 
point sources of EUV radiation (e CMa and [3 CMa) and low column densities 
are found. Both sets of models are constrained by LIC data obtained inside of 
the heliosphere, such as for pickup ions, He°, and anomalous cosmic rays. The 
PUI and ACR populations are seeded by neutral ISM flowing into the helio- 
sphere, and are subject to filtration losses in the heliosheath regions (Moebius 
et al., Zank et al., this volume). 

The ionization state of the pristine ISM outside of the heliosphere is ob- 
tained from the RT models, using slightly different constraints for the first 
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Figure 6.8. Variation of neutral densities between tlie Sun and cloud surface for Model 2 from 
Slavin and Frisch (2002). Shown are variations in neutral column densities between the Sun 
(Af(H°)=6.5x 10^^ cm-^) and cloud surface (7V(H°)=0) for H°, He°, Ne°, 0°, and N°. At the 
heliopause, n(He°)~0.015 cm~^, n(H°)~0.22 cm"^, n(e)~0.1 cm'^*. 



and second generations of models. The first generation of 25 models are con- 
strained by the local ISM towards e CMa, i.e. the sum of the LIC and blue 
shifted (BC) clouds, and by ISM byproducts inside the heliosphere. The ISM 
within ~l-2 pc towards e CMa is divided between two clouds, with projected 
HC velocities of 10 km s"^ (the BC) and 18 km s~^ (the LIC), and with a total 
A^(H°)< 10^^ cm^^. The LIC and BC are also observed towards Sirius (a 
CMa, 2.7 pc). Two models. Model 2 and Model 8, provide a good match to 
all data except for the cloud temperature inside of the heliosphere. The pre- 
dicted temperature at the heliosphere is ~3,000 K higher than found from in 
situ He° data, T= 6,400±340 K. A possible explanation for this difference is 
that the abundances of coolants such as C"*" are incorrect in the models, partic- 
ularly since A^(C"'") generally has large measurement uncertainties. A second 
mismatch occurs because predicted iV(Si"^^) values are lower than observed 
values, which suggests that interface models require additional processes such 
as a turbulent mixing layer. Models 2 and 8 predict that n(H°)~0.20 cm~^, 
n(e)~0.1 cm~^, and ionizations of x(H)~0.30, and x(He)'^0.49 are appropri- 
ate for the solar location (Slavin and Frisch, 2002, Frisch and Slavin, 2003). 
The mean cloud density to the downwind surface is (n(H°))Lic~0-17 cm~^. 
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Figure 6.9. Top: H neutral density vs. He neutral density, both at the Solar location for a large 
set of model calculations. The symbol colors indicate the assumed magnetic field strength in 
the cloud, the fill indicates the assumed HI column density, and the shape indicates the assumed 
temperature of the hot gas of the Local Bubble. Stars represent special parameter sets which 
do not fall on the grid of model parameters, but rather are chosen to better match the data. 
The density and uncertainties of He° are plotted as vertical lines Mobius et al., 2004. Note 
that a range of n(HI) values are consistent with the n(He°) measurements. The radiative transfer 
models that provide the best agreement with all available data on the ISM also yield a consistent 
estimate of n(H°)= 0.19 - 0.21 cm~^. Bottom: Same as left plot but for the H+ and He+ ion 
fractions. Note: Some of the differences not explained for by Af(H°), Tlb, or B are due to 
differences in the assumed total H densities. 



The second set of models is constrained by data on the LIC towards e CMa 
(excluding the BC) and in situ ISM. The best of these models are in good agree- 
ment with the in situ data, including cloud temperature data. However, these 
models require a high C abundance, possibly in conflict with solar abundances, 
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and fail to predict 7V(Mg+)/iV(Mg°) in the LIC. The second set of models re- 
quire that the LIC and BC have different ionization levels, which is similar to 
the findings of Gry and Jenkins (2001) but remains unexplained. LIC prop- 
erties predicted by these models, which are still under study, are 77,(H°)=0.19 
cm-3, x(H)=0.22, n(e)=0.06 cm-^, n(He°)=0.015 cm-^, and x(He)=0.37 

The results of the RT models that are significant for heliosphere studies are 
summarized in Figures 6.8 and 6.9. The ionizations of H, He, O, N, and Ne 
throughout the LIC are shown in Fig. 6.8, with variations of ~25% between 
the cloud surface and the Sun. Fig. 6.9 shows the extent of ionization states 
possible under equilibrium conditions in tenuous ISM. Minor variations in the 
physical assumptions input to the RT models result in a continuum of ioniza- 
tion levels that are in equilibrium for low density clouds like the LIC. De- 
pending on the RT model constraints, H ionization levels are 0.19-0.23, while 
He ionization levels are 0.32-0.53. These models show that ionization levels 
might vary between clouds in the CLIC, so that variations in the heliosphere 
boundary conditions are expected as the CLIC flows past the Sun (see §6.ij. 

6.5 Passages through Nearby Clouds 

The transition of the Sun from the near void of the Local Bubble, and into 
the stream of tenuous ISM flowing away from the Scorpius-Ophiuchus As- 
sociation, can be probed with the kinematics and column densities of nearby 
clouds, combined with models of the volume density n(H°). Specifically, cloud 
dimensions are assumed to be ocA^(H°)/(n(H°)). Radiative transfer models of 
tenuous ISM show that cloud ionization, x(H), varies with column density. 
Since the densities of individual cloudlets are not currently available, n(H°) 
determined by LIC radiative transfer models is extrapolated to other nearby 
clouds (§6.4). With this simple approach, we see that the Sun appears to have 
entered the CLIC within the past ~140,000 yeai^s. The Sun is located close to 
the downwind edge of the CLIC (Frisch, 1995). CLIC physical properties are 
discussed in §6.3. 

The most predictable passage of the Sun into an interstellar cloud is the 
epoch at which the Sun entered the LIC, the cloud now suiTounding the Sun. 
The LIC is the only cloud with an accurate 3D velocity vector, which is found 
from Ulysses measurements of interstellar He° inside of the solar system. 
The LIC HC velocity is 26.3±0.4 kms~^, and the downwind direction is 
A = 74.7 ± 0.5°, /3=-5.2±0.2° (ecliptic coordinates, Witte, 2004). Several 
estimates are given for the epoch of the Sun's entry into the LIC. Both the data 
and radiative transfer models, upon which the estimate relies, need further re- 
finement before answers are conclusive. We test several models of the poorly 
known LIC shape. The observed cloud velocity and the projected LIC velocity 
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towards the x Ori' 8.7 pc away and ~15° from the downwind direction, differ 
by ~2.4 km s~^, suggesting non-LIC gas is contributing (Table 6.4). 

The distance to the cloud edge is given by A^(H°)/(n(H°)), where (?i(H°)) 
is the average space density of H° in the cloud. The density of the LIC at 
the solar location is ?i(H°)~0.2 cm~^, and a mean density (n(H°))Lic~ 0.17 
cm~^ is found from radiative transfer models of the LIC in the downwind 
direction (§6.4). When the cloud structure is not fully resolved, or several 
velocity components are present, there may be gaps in the cloud that are not 
incorporated into the adopted value for (n(H°)). We present several estimates 
below for the entry of the Sun into the LIC, with results that differ by factors 
of three. However, the answers are consistent with the Sun's entry into the LIC 
sometime within the past ~47,000 years for (n(H°))=0. 17 cm"'^, and possibly 
quite recently within the past thousand years. 

The difference between the LIC and upwind gas velocities, of <, 2 km s~^, 
including towards the nearest star a Cen, suggest that the properties of the 
cloud surrounding the Sun will change rather soon, in ^^3,800 years. 

6.5.1 First Encounter of Sun with the Local ISM 

The projected LIC velocity and A^(H°) towards nearby stars in the down- 
wind direction (Table 6.4), indicates that the Sun entered the CLIC sometime 
within the past 44,000-140,000/?ii^/,o.2 years, where n///,o.2 is the average 
cloud density in units of 0.2 cm~'^ (§6.4, Fig. 6.9). The value n//7^o.2~l is 
reasonable for nearby tenuous ISM, and is consistent with radiative transfer 
models of the LIC. This ti'ansition of the Sun out of the deep vacuum of the 
Local Bubble, and into the higher density CLIC ISM, would have been accom- 
panied by the appearance of neutral interstellar gas, pickup ions, anomalous 
cosmic rays, and dust in the heliosphere. The geological record of cosmic ray 
radioisotopes should have sampled this transition (see the discussions of cos- 
mic rays and the spallation product radioisotopes in Chapters 9, 10, and 12). 

The simplest estimate for the date of entry of the Sun into the CLIC is to 
ignore non-radial motions for all clouds and look only at the nearby stars in 
the anti-apex direction. ISM towards the three stars a CMaA, a CMaB, and 
e CMa are useful for this estimate, giving A^(H°)=5-7 x 10^^ cm^^ (based on 
data and models in Hebrard et al., 1999, Frisch and Slavin, 2003). The distance 
to the cloud edge in the anti-apex direction is (xnHifl.2~^ pc The most distant 
of the two cloudlets seen in these sightlines is a blue-shifted cloud, with a 
heliocentric radial velocity of 9.2-13.7 km s^^ (Table 6.1). We then find that 
the Sun entered the CLIC within the past (59,000-120,000)/ni//,o.2 years. 

Using instead the ISM towards the stars Ori ^rid a Aur, which are near 
the heliocentric downwind direction, and for a distance to the cloud edge given 
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by N(}l°)lnHifl.2, the observed radial velocities (Table 6.4) indicate that the 
Sun entered the CLIC (44,000-140,000)/ni^/_o.2 years ago. 

For these estimates, possible non-radial cloud motions and gaps between 
clouds are ignored. The cloud gas filling factor in the downwind direction is 
fr^O.26 in §6.3. If the downwind CLIC consists of a series of wispy clouds 
with similar velocities, then nni, o.2'^0.2 cm""^, and the first solar encounter 
with the CLIC could have occurred up to 4 x 10^ years ago based on the 
Ori data. 

6.5.2 Entering the Local Cloud 

In principle the entry date of the Sun into the surrounding cloud, LIC, can 
be determined precisely since the full 3D velocity of the LIC is known from 
observations of He° inside of the heliosphere (Table 6.1). The LIC shape is 
uncertain because components at the LIC velocity are not resolved in many 
sightlines. Although there may be high density contrasts from low speed waves 
in the LIC, for instance acoustic waves will have velocities ~7 km s~^, we take 
the approach that local clouds are defined by the component structure and look 
at several geometric LIC models to estimate the epoch the Sun entered the LIC. 
ISM data for downwind stars are summarized in Table 6.4. 

The simplest estimate (A) uses A^(H°) and the LIC velocity towards 
Ori, ^^15° from the downwind LIC direction. Estimate (B) is based on the 
closest observed star in the downwind direction, a CMa, which is ~45° from 
the downwind direction, combined with the assumption that the normal to the 
downwind cloud surface is parallel to the LSR LIC velocity. This assumption 
results in an entry epoch that varies with the assumed LSR (§6.1.2), and ig- 
nores data for more distant downwind stars. Estimate (C) approximates the 
downwind LIC surface as a flat plane, defined by the distance of the LIC edge 
towards any three stars from Table 6.4, moving through space at the LIC ve- 
locity. The final estimate (D) relies on the Colorado LIC model (RLO), which 
however does not incorporate the most recent data for the LIC, but can be ex- 
pected to substantially improve in the future. 

A. The values for the LIC A'^(H°) and the projected LIC velocity towards x^ 
Ori (Table 6.4) indicate that the Sun entered the LIC ^^40, 000/ uhi, 0.2 years 
ago, where nHifl.2 is the mean n(H°) to cloud surface in units of 0.2 cm^'^. 
Radiative transfer models indicate n^^/ o.2~0.85 (§6.4). The observed cloud 
velocity towards Ori differs from the projected LIC velocity by 2.4ibl.O 
km s~^, so that additional cloud gas is present in this sightline that is not at the 
LIC velocity. In addition, the turbulent contribution to the line broadening is 
i = 2.38;[J;}^ km s-i towards x^ Ori (§6.1.3, Redfield and Linsky, 2004b), 
indicating unresolved cloudlets may be present. 
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B. For this estimate we use LIC data towards the nearest star in the down- 
wind direction for which ISM data are available, a CMaAB (Sirius). The LIC 
downwind surface is assumed to be oriented such that the surface normal is 
parallel to the LIC LSR velocity, based on the model of Frisch (1994, Fig- 
ure 6.10). Note that the "downwind" nomenclature is referenced with respect 
to the Sun, and the "downwind surface" is the leading edge of the LIC as it 
moves through the LSR. Since the LIC surface is referenced to the LSR, the 
solar apex motion is a variable (Table 6.1, §6.1.2). The Sun then encountered 
the LIC within the past ~13,500/n/f/ 0.2 years or ~ lO'^/rij^/ 0.2 years, for the 
Hipparcos or Standard solar apex motions, respectively. The present distance 
to the LIC downwind surface is ~ 0.01 — 0.3 pc for this model. The velocity 
of most of the interstellar gas observed towards x Ori ^rid a Aur differs by 
a small amount, 0.7-1.4 km s~^, from the projected LIC velocity (Table 6.4), 
and therefore is located in a separate cloud. Because the LIC also extends less 
than ~ 10"^ AU in the upwind direction, this model implies that the LIC is fil- 
amentary or sheet-like, similar to global low density ISM (§6.6). This model 
results in an interstellar- magnetic field direction, i?is, which is approximately 
parallel to the cloud surface (perpendicular to the surface normal). 

Figure 6.10. SEE fig6.pdf. The distribution of nearby ISM is shown based on the assumption 
that the LSR motion of the LIC is parallel to the normal of the cloud surface (model B, §6.5.2). 
The Sun and LIC motions are in the LSR. Figure from Frisch (2000), courtesy of American 
Scientist. 



Table 6.4. Stars Sampling the Downwind LIC 



Object 


l,b 


Dist. 


^observed 




logiV(H°) 




(°,°) 


(pc) 


(km s~i) 


(km s-i) 


(cm^'^) 


a Aur 


162.6, 4.6 


12.9 


21.5±0.5 


23.0±0.3 


18.26 


LIC He° 


183.3,-15.9 





26.3±0.3 


26.3±0.3 




Ori 


188.5, -2.7 


8.7 


23.1 ±0.7 


25.5±0.3 


17.80 


aCMi 


213.7, 13.0 


3.5 


20.5±1.0, 24.0 


19.6±0.2 


17.90 


a CMaA 


227.2, -8.9 


2.7 


13.7, 19.5±0.3 


19.1±0.2 


17.25 


aCMaB 


227.2, -8.9 


2.7 


11.7, 17.6±1.5 


19.1±0.2 


17.63 


eCMa 


239.8, -11.3 


132 


9.2, 16.2±1.5 


15.1±0.2 


~17.54 



Notes: The LIC velocity component is in italics. V^jq is the LIC velocity 
projected towards the star. All velocities are heliocentric. The column densities 
and data are from Witte (2004), Hebrard et al. (1999), Frisch & Welty (2005), 
RL, GJ, SF02, andFS. 



C. The third estimate models the downwind cloud surface as a locally flat 
surface moving with the LIC velocity, based on the model in Mueller et al. 
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(2005). We then use any three stars in Table 6.4 to define this plane, and 
the time at which the Sun crosses this plane marks the entry of the Sun into 
the Lie. For this model, the Sun entered the LIC between 33,000/n/// 0.2 to 
36,000/nH/,o.2 years ago. 

D. The final estimate relies on the Colorado LIC model, based on density 
(n(H°))=0.1 cm~^ and a LIC velocity determined from UV absorption lines 
towards nearby downwind stars (RL). This UV LIC vector differs by ^^0.6 
km s~^ in velocity, and ^^0.8° in direction, from the LIC velocity given by 
data on interstellar He° inside of the solar system (Witte, 2004). The Colorado 
Model Web calculator gives a distance to the downwind surface of 4.5 pc, in- 
dicating the Sun entered the LIC 170,000 years ago for a relative Sun-LIC 
velocity of 26.3 km s~^, or alternatively ~85, 500/^/^/ 0.2 years ago. The col- 
umn density towards the downwind direction is predicted to be A^(H°)=10^^'^^ 
cm~^, however, which differs substantially from recent values for Ori (Ta- 
ble 6.4). We therefore disregard this estimate, but expect further improvements 
in the Colorado model to significantly increase our understanding of the LIC 
morphology. 

6.5.3 Future Cloud Encounters 

The best limits on the distance to the upwind edge of the LIC are found 
from observations of 36 Oph, 6 pc away and ~ 16° from the heliosphere nose. 
The ISM velocity towards 36 Oph is -28.5±0.6 km s~^, which corresponds 
to the G-cloud velocity (Table 6.3). The limit on A^(H°) at the projected LIC 
velocity of 25.7 km s^^ is iV(H°)< 6 x 10^^ cm^^ (Wood et al., 2000b). 
These numbers indicate that the Sun will exit the LIC sometime within the 
next ~3700/n/// 0.2 years, and that the distance to the upwind LIC surface is 
<0. 10 pc. Support for this result is provided by the ISM velocity towards the 
nearest star a Cen, which is ^^50° from the HC nose direction. Towards a 
Cen, the observed cloud velocity and the projected LIC velocity differ by '^l 
km s~^ (Landsman et al., 1984, Linsky and Wood, 1996). Column densities of 
A^(H°)= 17.6-1 8.0 cm~^ were originally found towards a CenAB (Linsky and 
Wood, 1996), which indicate that (n(H°)) may be similar to the LIC density of 
^^0.2 cm^^. Later interpretations of the data favor the low end of this column 
density range (Wood et al., 2005), indicating that the sightline is only partially 
filled with warm gas. 

Possibilities for the next cloud to be encountered by the Sun include the G 
and Apex clouds (Table 6.3, Frisch, 2003). Both clouds have LSR upwind 
directions, {i, 6)=(348ib3°, 5±5°), that are approximately perpendicular to the 
LSR solar apex motion (§6.1.2). Observations of a Cen and a Aql show that 
the G-cloud is within 1.3 pc, and the Apex cloud is within 5 pc of the Sun, 
giving upper limits for an encounter date of ^^45,000 and ~ 175,000 yeai^s. 
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respectively. The column densities of the G-cloud are ~3 times larger than 
Apex cloud column densities, which are A^(H°)~ lO^'''^'' cm~^ towai^ds 70 
Oph and a Aql. Typical G-cloud column densities are A^(H°)~ 10^^'^^ cm~^. 

Both clouds are inhomogeneous. Only the Apex cloud is seen towards AU 
Mic (d=9.9 pc, ^=12.7°, 6=-36.8°). Only the G-cloud is seen towards o Ser 
(d=13 pc, £=13.3°, 6=9.2°) and v Ser (d=59 pc, £=10.6°, 6=13.5°). Both 
clouds are seen towards a Aql (d=5 pc, £=47.7°, 6=-8.9°) and a Oph (d=14 
pc, £=35.9°, 6=22.6°). 

The G-cloud Ca+ component (iV(Ca+)=1.58x 10^^ cm"^) towards a Oph 
is extraordinarily strong for a nearby star, and it implies a density of n(H°)>5 
cm~^ in the G-cloud, which is 15 times the density of the LIC (Frisch, 2003). 
This estimate is found by using the G-cloud limit of iV(Ca+)/iV(H°)< 7.1 x 
10"^ towards a Cen in order to estimate A^(H°) for the G-cloud in the a Oph 
sightline, and by assuming the G-cloud is entirely foreground to a Cen, giving 
a distance limit of 1.3 pc. However, should the G-cloud temperature be colder 
towards a Oph than towards a Cen, or should a significant amount of the G- 
cloud be beyond the distances of 36 Oph and a Cen, then the requirement for 
a high density becomes diminished. 

6.6 The Solar Environment and Global ISM 

The global properties of ISM in the solar neighborhood offer hints about 
potential variations in the ISM that might impact the heliosphere. 

Over the past ^^50 years, radio and UV data have slowly disproved simple 
ISM models consisting of gravitationally clumped clouds embedded in a sub- 
strate of warm or hot low density gas. Evidence supporting a filamentary nature 
for low column density ISM, with iV(H)< 10^^ cm~^, is provided by 11° 21 
cm observations of high latitude H° filaments, and also by UV observations 
that are interpreted as high space density (n(H°)>10 cm~'^), low column den- 
sity ( iV(H°)< lO^'^ cm"2 and iV(H+)< 10^^ cm"^) nearby gas (Hartmann 
and Burton, 1997, Welty et al., 1999, Welty et al., 2002). Towards 23 Ori, 
warm low velocity diffuse gas has densities in the range n(H°)= 15-20 cm~^, 
and cloud thicknesses, L, are 0.7-0.9 pc. The densities of cold clouds, with 
temperatures <100 K, are slightly lower, n(H°)=10-15 cm~^. Diffuse ionized 
gas is present at intermediate velocities, |Vlsr| > 20 km s~^, n(e)<5 cm~^, 
and L= 0.0006-0.04 pc. Thin diffuse ionized clouds, with n(e)~0.17 cm~^, 
L= 0.2-2.7 pc, and similar velocities, are observed towards Orion stars spread 
over >15°, which indicates that these clouds must be thin sheets or filaments. 

Winds of massive stars and supernova in young stellar associations eject 
energy into the surrounding ISM, and evacuate cavities in the ISM that are 
suiTounded by irregular shells of WNM. These shells, which may overlap or 
be filled with X-ray plasma, are traced by H° 21 -cm emission, and explain 
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some of the ISM filamentary structure (MacLow and McCray, 1988, Cox and 
Smith, 1974, Dickey, 2004, Heiles and Troland, 2003b). The CLIC appeal's to 
be a fragment of a superbubble shell from one epoch of star formation in the 
SCA, and a Oph may be in a direction that is tangential to this shell (Frisch, 
1995). 

6.6.1 Warm Ionized and Partially Ionized Gas 

Interstellar plasma and neutrals interact with the heliosphere in fundamen- 
tally different ways, and small variations in the interstellar radiation field con- 
vert neutral gas to plasma. Nearby WPIM was discussed in §6.3.1, and we 
now compare warm interstellar plasma over 100 pc scales and show that the 
CLIC and global observed diffuse ionized gas are similar. We conclude that 
dense ionized ISM has not been, and will not be, part of the immediate solar 
environment over short timescales. However, diffuse WPIM similar to the LIC 
and CLIC is much more widespread. 

Dense fully ionized ISM is not found close to the Sun, although diffuse 
ionized gas surrounds many hot stars at distances of >150 pc, such as j3 CMa 
bordering the giant ionized Gum Nebula, j3 Cen, which energizes the interior 
of Loop I, and several stars in Upper Scorpius. Generally H II regions surround 
hot O, early B, and white dwarf stars that emit strongly in the EUV, A <912 
A. The solar path is unlikely to traverse an H II region surrounding an 0-B 1 
star over timescales of ±4 Myrs because there are no hot stars or dense H II 
regions nearby. The closest H II region surrounds the high latitude star a Vir, 
80 pc away. White dwarf stars are relatively frequent near the Sun. Over 25 
white dwarf stars with surface temperatures T > 10^ K are found within 20 pc, 
and > 40% of these stars are hotter than 15,000 K. Small Stromgren spheres 
will surround the hottest white dwarf stars, but the densities will be very low 
for nearby stars (Tat and Terzian, 1999). 

Diffuse ionized gas is widespread, in contrast to dense ionized gas. It is 
traced by weak optical recombination lines such as the Ha line, N II, He I, 
He II, O I, and Si II (Haffner et al., 2003, Reynolds, 2004), UV absorption lines 
(Welty et al., 1999, Holberg et al, 1999, Lehner et al., 2003, GJ), and EUV ob- 
servations of He° and He+ towards white dwarf stars (§6.3). Pulsar dispersion 
measures trace the interactions between pulsar wave packets and electrons (e.g. 
Taylor and Cordes, 1993, Cordes and Lazio, 2002, Armstrong et al., 1995). 
These data reveal a low density, ionized component, n(e)~0. 1 cm~^, which 
fills ~20% of the volume of a 2-kpc-thick layer around the Galactic disk. The 
diffuse H"*" has been mapped in the northern hemisphere by WHAM, however 
the LIC emission measure is an order of magnitude below the WHAM sensi- 
tivity (e.g. Reynolds, 2004). About 30% of diffuse H° and H+ gas are found 
in spatially and kinematically associated clouds with densities ^■0.2-0.3 cm~^ 
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(Reynolds et al., 1995). If the H° and H+ were spatially coincident in these 
clouds then x(H) would be ~ 0.4, but Reynolds et al. argue that the H° and 
are spatially sepai^ated. Despite the evidence from spatial maps and line widths, 
however, we believe that the existence of some true WPIM (perhaps at lower 
ionization levels) cannot be ruled out in these regions. These clouds ai^e associ- 
ated with large filamentary or sheet-like structures, and exhibit cloud densities 
and hydrogen ionization levels that are similar to LIC values. In other direc- 
tions. Ha line emission indicates that H is highly ionized in the WIM, x(H)~l. 
Emission at 588 A from He° shows that for the WIM x(He)~0.3-0.6. Ratios 
of x(He)/x(H) are ~0.3-0.6 are found, versus ~1.6 for the LIC, implying that 
the LIC radiation field is harder than for most of the WIM, perhaps because 
of lower column densities. The WIM has a continuum of temperatures, 5000- 
10,000 K, and ionization levels, as demonstrated by emission in [Nil], [Ol], 
and [S II] lines that trace temperature and ionization. 

The interstellar plasma is structured over all spatial scales, and these plasma 
clumps appear to be similar- to the LIC. The plasma is partially opaque to ra- 
dio emission in the energy range 0.1 to 10 MHz, and as a result synchr^otron 
emission probes the WIM dumpiness (Kulkarni and Heiles, 1988, Peterson 
and Webber, 2002). Models for the propagation of low frequency synchrotron 
emission through a clumpy WIM fit the radio data better than do propagation 
models, that assume a uniformly distributed parallel slab model for the WIM. 
These clumps have density n(e)~0.2 cm~'^, fill 8%-15% of the disk, and have 
a free-free opacity at 10 MHz that is consistent with WIM temperatures as low 
as 4500 K. The dispersion, refraction and scintillation of pulsar wave packets 
indicates that the WIM is turbulent over scale sizes 10-2 - 102 AU, or even 
larger, ~ 10^ pc, if Faraday rotation and plasma density gradients are included 
(Armstrong et al., 1995). 

The WIM properties are consistent with predictions of the radiative transfer 
models, which find x(H+)=0. 15-0.35, T ~ 5, 000 - 9, 000 K, and n(H°)=0.2- 
0.3 cm^'^ for low column density clouds (Fig. 6.9, Slavin and Frisch, 2002). 
For very low column densities, <10^^ cm^^, a continuum of ionization and 
temperature levels are expected from variations in the cloud column density 
and the EUV radiation field. The WNM and WIM in the CLIC may show 
similar variations, if the A^(H°)~10^^ cm-^ clouds represent clusters of lower 
column density objects. The uncertain role of turbulence in cloud evolution 
allow this possibility. 

6.6.2 Neutral Interstellar Gas and Turbulence 

Radio and UV observations of low A^(H°) cloudlets in the solar vicinity in- 
dicate that dense, low column density, cloudlets may be hidden in the CLIC 



The Solar Environment and Global ISM 



45 



flow. Radio 21 -cm data indicate ISM of this type is widespread, and has col- 
umn densities and kinematics similar to CLIC ISM. 

Two populations of interstellar clouds, WNM and CNM, are found based 
on observations of the collisionally populated H I hyperfine 21 -cm line towards 
radio-continuum sources. The recent Arecibo Millennium Survey used an ON- 
OFF observing strategy to survey the properties of H° clouds towards ~80 
radio continuum sources, including emission and and opacity profiles, which 
were then fit with ~375 Gaussian components. The detailed attention paid to 
sources of noise and line contamination yielded a data set that provides a new 
perspective on the statistics of low A^(H°) CNM and WNM clouds (Heiles and 
Troland, 2003a, Heiles and Troland, 2003b). 

The CNM and WNM components observed in the Arecibo survey show that 
^^25% of the ISM H° mass is contained in clouds traveling with velocities 
> 10 km s~^ through the local standard of rest. When the solar apex motion of 
13-20 km s~^ is included, this means that Sun-cloud encounters with relative 
velocities lai^ger than 25 km s~^ are quite likely, as we might have guessed 
because of the LIC heliocentric velocity of -26.3 km s~^. 

For 21 -cm absorption components with spin temperature Tg, the opacity is 
r ocA^(H°)/Ts and the line width is FWHMocfjv, where is the dispersion in 
cloud velocity. Therefore the turbulent (A14urb) and thermal contributions to 
(7v can be distinguished for the CNM. Spin temperatures, Tg, for the CNM are 
typically ~20-100 K. 

The Arecibo survey discovered that 60% of the diffuse ISM is WNM, and 
~40% is CNM. Median column densities for the WNM and CNM towards 
sources with \h\ >10° are iV(H°)=1.3 x 10^^ cm-^ and 5.2 x 10^^^ cm-^ 
respectively. The finding that CNM column densities are typically lower than 
WNM values was unexpected. The column density weighted median spin tem- 
perature for the CNM is 70 K, although minimum temperatures are 20 K or 
lower. Upper limits on WNM kinetic temperatures are 15,000 K or higher, 
with a typical value of 4,000 K. The turbulent velocities derived for CNM 
components indicate that the turbulence must be highly supersonic, with a tur- 
bulent Mach number 1/tuib = "^^^tuThl ^^3, where is the sound speed 
and the factor of 3 converts to thixe-dimensional turbulence (Heiles, 2004b). 
The densities of the CNM observed towards 23 Ori are > 10 cm~^ (Welty et al., 
1999). 

The velocities of CLIC, WNM, and CNM components are compared in Fig. 
6.11. Only high latitude WNM and CNM data with |6| >25° are included, 
in order to avoid distant gas near the galactic plane. Some of the CNM and 
WHM components are found at velocities lower than -25 km s~^, which may 
be partly due to infalling H° gas at high latitudes. Otherwise, the CLIC kine- 
matics and column densities overlap generic WNM and CNM values. 
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Figure 6.11. SEE figures figlla.pdf and figllb.pdf. Left: CLIC velocities, for d<50 pc, are 
plotted against CNM and WNM velocities for sightlines with |6| >25° in the Arecibo Millen- 
nium Survey Heiles and Troland, 2003b. The Arecibo clouds are within -^425 pc if the H° layer 
thickness is 185 pc. The CLIC iV(H°) is based on D°/H°= L5 x 10~^ Ca+/H°= 8.3 x 10"^ 
or Mg+/H°= 6.7 x 10~®, and LSRgtd is used. Right: This composite profile shows that some 
CNM components can be lost in the WNM at UV resolutions. The solid line shows a D° 1215 
A absorption component formed by two clouds at the same velocity (0 km s^^) and column 
density {N(D°)= 10^^ " cm^^), but different temperatures (50 K and 7,000 K). The dashed 
line shows WNM only, with iV(D°)= lO" " cm"^ and T = 7000 K. This theoretical profile 
incorporates a Gaussian instrument resolution of FWHM 3 km s~^ appropriate for UV data. 
Turbulent broadening is not included. 



Of special interest are very low column density cold CNM, or "tiny" clouds, 
found for a broad range of velocities. An example is a component with Ts=16.6 
K, iV(H°)=2 X 10^^ cm~2 and velocity -40.6 km s~^ observed towards 3C 
225. Another example is a component with Ts=43 K, iV(H°)= 3 x 10^^ 
cm~^, and velocity -4.6 km s~^ observed towards 4C 32.44. In addition, CNM 
components with iV(H°)~10^^ cm"^, Ts~30-90 K, and velocities -30 to 
km s~^, are confirmed towards several other sources (Stanimirovic and Heiles, 
2005, Braun and Kanekar, 2005). 

CNM ISM is detected in the UV towards 23 Ori, with low column densities, 
weak ionization, and low temperatures: A^(H) ^ 2 x 10^^ cm~^, rs~60-150 
K, and ~1% ionization. Cold clouds with lower column densities may be 
difficult to resolve at UV resolutions. The densities of these clouds are found 
from pressure considerations using C I and C II fine structure lines (Jenkins, 
2002), and yield nH'^10-15 cm^^ and a cloud thickness of ~0.5 pc (Welty 
et al., 1999). These components are formed in the Orion-Eridanus Bubble 
shell. Similar features are observed in other Orion stars. This CNM is either 
filamentary or sheet-like, and it is a plausible template for the CNM in Fig. 
6.11. 

The evident crowding of absorption components in velocity space may in- 
troduce problems with the interpretation of UV data (§6. 1.3). CNM and WNM 
clouds with low A^(H°) and similar velocities would be unresolved at a nomi- 
nal 3 km s~^ UV resolution, as shown by the theoretical composite D° profile 
plotted in Fig. 6.11. This profile is synthesized from the expected profiles of 
two clouds at km s~^, temperatures 7,000 K and 50 K clouds, and each with 
A^(H°)=10^^ cm~^. An instrumental resolution of 3 km s~^ is included, but 
turbulence is not. Resolution problems are worse for heavy elements, since the 
line FWHMoc mass~^^^. Therefore local low A^(H°) cold cloudlets can not be 
ruled out, and if present would be traversed by the heliosphere in ~ 1-20,000 
years. 
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6.6.3 Magnetic Fields 

Interstellar magnetic fields interact strongly with the heliosphere, and mod- 
ify the propagation of galactic cosmic rays in space. Unfortunately the uni- 
formity and strength of the interstellar magnetic field (Sis) at the Sun are un- 
known. We showed that magnetic field strengths of i?is~0.6-3.7 i^lG would 
satisfy several assumptions about the equilibrium of local ISM (§6.3.4). How- 
ever, the only direct data on the nearby magnetic field are the starlight polar- 
ization data attributed to weakly aligned interstellar grains in the nearest ^■SS 
pc (Tinbergen, 1982, Frisch, 1990). More recently, it appears as if these grains 
are instead trapped in the magnetic wall of the heliosheath (Frisch, 2005). The 
global nearby Bi^ thus must be used to obtain insight on i?is near the Sun. 

The properties of the interstellar magnetic field, including the weak spa- 
tially uniform {Ba) and stronger random {B^) components, are known for both 
CNM and WIM (e.g. Beck, 2001, Crutcher et al., 2003, Heiles, 2004a, Heiles 
and Crutcher, 2005). The CNM median field strength of 6.0±1.8 ^jlG is de- 
termined from data on the Zeeman splitting of the H I 21 cm line (Heiles and 
Troland, 2005), but may not apply to the CLIC where CNM is not yet de- 
tected. Although i?is~ n for densities larger than 10^ cm~'^, at lower densities 
i?is strengths are uncon^elated with cloud density, in violation of simple flux 
freezing assumptions. 

Pulsar data, synchrotron emission, and starlight polarization data indicate 
that the ratios of the random {B,-) and uniform {B^) components of the lo- 
cal field are BJBj-^O.b-l.O. Synchrotron emission of external galaxies indi- 
cate that the ratio BJB^ increases in interaim regions such as suiTounding the 
Sun, but decreases in spiral arm regions as the magnetic fields become more 
random. Comparisons of magnetic field strengths determined from polarized 
synchrotron emission and pulsar dispersion measure observations, which in- 
teract with Bis and the WIM, give the total local field strength and the rela- 
tive strengths of i?u and Bj-. Synchrotron emission gives a total field strength 
locally of i?is~ 6 ± 2 /iG, and the uniform component strength -Bu~ 4 fiG. 
Faraday rotation from the propagation of the oppositely polarized circular com- 
ponents produces a polarization phase shift, given by the rotation measure RM 
= i?l|(x)ne(x)dx, which responds to the parallel component of the magnetic 
field, weighted by the electron density, nc{x), and integrated over path- 
length L. Bis is directed away from the observer for RM<0. Wave packet 
dispersion in the WIM gives the dispersion measure DM = nc{x)dx. The 
mean local field, weighted by n(e), is found from the ratio RM/DM and is 
6.5 fiG, where B^^B\\. 

Starlight polarization data show that the local Bis is relatively uniform, and 
not strongly strongly distorted or tangled by turbulence. Locally the total field. 
Bis, follows the pitch angle, 9.4°, of the local spiral arm, and has a curvature 
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radius of ~7.8 kpc centered on ^=344.6°. The starlight polarization vectors 
converge near ^~80° and ~270°, where the random component of the field, 
Bt-, results in a polarization angle that is rotated over all angles. 

A large excursion from is seen towards within 150 pc, towards Loop 
I, where both polarized starlight and synchrotron emission indicate a mag- 
netic field conforms to the superbubble shell around the Scorpius-Centaurus 
Association. Surveys of B\\ in morphologically distinct objects such as the 
Loop L Eridanus, and Orion supershells, show that the magnetic field is strong, 
-By > 5 ^G, and ordered throughout these shell features. 

6.6.4 Origin of Local ISM 

Several origins have been suggested for the CLIC. CLIC kinematics and 
abundance patterns led to the suggestion that the CLIC is part of a superbub- 
ble shell surrounding the Scoipius-Centaurus Association (Frisch, 1981, Frisch 
and York, 1986, de Geus, 1992, Frisch, 1995). The SCA is surrounded by a 
well-known supershell formed by epochs of star formation, winds from mas- 
sive stars and supernovae (Weaver, 1979). Using superbubble expansion mod- 
els and the star formation rates in the SCA, Frisch modeled the CLIC as a 
fragment of the shell from the formation of the Upper Scorpius subgroup (~ 4 
Myrs old) of the SCA, which has expanded away from the high-density ISM 
close to the association to the lower pressure regions around the Sun (Frisch, 
1995). A related suggestion attributes the origin of the CLIC to Rayleigh- 
Taylor instabilities in the shell of the younger Loop I supernova remnant (~0.5 
Myrs old), which has propelled ISM towards the Sun (Breitschwerdt et al., 
2000). A third origin suggests that the CLIC originates as a magnetic flux tube 
that is originally captured in the inner wall of the Local Bubble suiTounding 
the Sun, and is driven back into the Local Bubble interior by magnetic tension 
(Cox and Helenius, 2003). This scenario requires that the CLIC velocity and 
magnetic fields are parallel, which appears unlikely because of the asymmetric 
capture of small interstellar- dust grains in the heliosheath (Frisch, 2005). Each 
of these suggestions is consistent with the CLIC kinematics, but should have 
different consequences for short-term variations in the solar environment. 

6.7 Summary 

In this chapter we have tried to present the essential information required to 
understand short-term variations in the galactic environment of the Sun, and at 
the same time make an educated guess about the times that the Sun has tran- 
sitioned between different types of ISM for the past 10'^ — 10^ years. It is our 
hope some day that the geological radio isotope records can be searched for 
the signatures of these transitions, although earlier searches have been incon- 
clusive (e.g., Frisch, 1981, Sonett et al., 1987, Frisch, 1997). 
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Both the Sun and ISM move through our local neighborhood of the Milky 
Way Galaxy, which is described by the LSR rest velocity frame (with its own 
uncertainties). The Sun has been in the nearly empty space of the Local Bub- 
ble interior for several million years. Sometime in the past ~140,000/nH/,o.2 
years, the Sun exited the nearly empty Local Bubble and entered an outflow 
of ISM with an upwind direction towards the Scorpius-Centaurus Association. 
Even more recently, ~40,000/n/f/ 0.2 years ago, the Sun entered the cloud now 
surrounding us, the LIC. Within ~4,000 years the Sun will exit the LIC. What's 
next? The flow of local ISM past the Sun is bearing down on the Sun at a rel- 
ative velocity of ~28 km s^^, the cloudlets in this flow are expected to pass 
over the Sun more often than roughly once per ~70,000 years. The next cloud 
the Sun is most likely to traverse is either the G-cloud or the Apex Cloud. This 
cloud must be either either inhomogeneous or denser by a factor of ~30 than 
the LIC. 

Most of this neai^by ISM is warm partially ionized gas, not dissimilar to 
warm ISM observed elsewhere, except that column densities tend to be lower. 
The properties of this gas are mainly consistent with predictions of radiative 
transfer models, which then suggests that the local gas is close to equilibrium. 
If viewed from elsewhere in the galaxy, the ensemble of ISM close to the Sun 
could easily pass as warm neutral material such as observed by Hi 21 cm 
data. The general characteristics of nearby ISM are similar to the global warm 
neutral material, except that column densities are very low so that H° ionizing 
photons penetrate to the cloud cores. The physical properties of the CLIC 
vary both between and within clouds. The low column densities of this ISM 
lead to ionization variations inside of a cloud. As the Sun passes between 
clouds, the interstellar ionization, dynamic pressure or the thermal pressure 
may vary, possibly all at once. As a consequence, the boundary conditions of 
the heliosphere may experience dramatic variations as the Sun moves thi^ough 
the CLIC. The efforts to understand these variations have just begun. 

One of the more intriguing possibilities is that tiny cool ISM cloudlets ai^e 
undetected in the nearby warm partially ionized gas (WPIM). Such structures, 
with A^(H°)~ 10^^ cm~^, are seen elsewhere, although infrequently. A new 
class of cloud models, perhaps based on macro-turbulence, will be required if 
such structures exist locally. 

Most of our conclusions have substantial inherent uncertainties. The most 
prominent is that all distances, such as those to cloud edges, are based on 
A^(H°), which is moderately well known for most sightlines discussed here, 
and an assumption about the mean H° space density, (?i(H°)). We have a very 
good idea of the value for (n(H°)) in the LIC because of extensive develop- 
ment of radiative transfer models for the LIC. Our confidence in the models 
follows from their ability to predict both ISM densities inside of the helio- 
sphere, including the Ulysses satellite values for n(He°), and ISM column 
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densities towards nearby stars. Since column densities are too small locally 
for measurements of the fine-structure lines that yield n(H°), we have assumed 
that the LIC value for (n(H°)) applies elsewhere in the CLIC. In addition, the 
transition epochs for the Sun are calculated by assuming that there are no gaps 
between individual cloudlets, as defined by velocity. The LIC (n(H°)) value is 
probably representative of the downwind CLIC, because the Sun is so close to 
the downwind edge. However, we can not rule out gaps between the cloudlets, 
and such gaps would alter estimates of the epoch the Sun entered the CLIC. 

One of the most significant improvements will be to make good maps of 
nearby ISM. Such maps would require a new generation of space instrumenta- 
tion, capable of high resolution, R~500,000, and high signal-to-noise observa- 
tions of both bright and faint objects over the spectral range 912-3,000 A. We 
also need more precise measurements of the diffuse interstellar radiation field, 
throughout the full spectral interval of 2000 A to 0.1 keV. Such data would 
provide a strong constraint on the interface emission and reduce uncertainties 
in the Mg° photoionization rate. 

In this chapter we have reviewed the short-term variations in the galactic 
environment of the Sun. Our galactic environment affects the heliosphere, 
and by analogy the astrospheres of neai^by stars, so that the space trajectory 
of a star is a filter for other planetary systems with conditions conducive to 
stable climates for exoplanets (Frisch, 1993). As other papers in this volume 
demonstrate, the ISM-modulated heliosphere has a pronounced effect on the 
cosmic ray flux in the inner heliosphere, which in turn appeai^s to affect the 
climate. Similar- processes will affect the climates of planets outside of the 
solar system. Such possibilities make the study of short-term variations in the 
galactic envir-onment of the Sun highly topical. 
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